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Least  Squares  Kit  of  a Nonlinear  Curve 

The  I o l loving  discussion  concerns  a method  ior  I 1 tiding  the  value.  of 
M parameters  ' , . * , in  the  formula 
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which  render  the  formula  a 'best  lit"  to  N measured  points,  having 

ordinates  Y|,  Y j , Y , '*’  Y and  abscissae  X , X , X , *“X  . The 
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formula  will  he  considered  a best  fit  If  t lie  goodness  of  fit  parameter 
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is  an  absolute  minimum  with  respect  to  the  M parameters  {>•.  The  parameters 

Wj  are  weighting  factors,  reflecting  measurement  reliability,  and  in 

( 

tii is  case  all  are  unity. 

The  necessary  criterion  for  A to  be  a minimum  with  respect  to  the  M 
parameters  {\}  is  that 
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When  the  formula  for  F is  not  linear  in  the  M parameters  1 A > , Equations(3) 
will  not  be  linear  in  the  same  M unknown  parameters.  However,  a set  of 

linear  equations  can  bt  obtained  from  Equations  (3)  by  approximating  eacli 

| • 

partial  derivative  ^ by  the  linear  terms  of  a Taylor  series  expansion 

m 

about  an  initial  point  {>o). 
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Tin  partial  derivative  terns  in  I.  ;uat  ion  (4)  i .in  lie  expressed  ■ • t i •• 


Equations  (4),  (5)  and  (h)  can  now  l>e  combined  to  yield 

* i:  + A (A  - •,  ) 0 (m  * 1,  2,  3,  ’•!) 
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where 
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The  column  matrix  of  improved  values  of  the  M parameters  (•  can  he 
calculated  by  inverting  Equation  (8). 

{X}  = fxo}  - a"1  { C } (11) 


The  process  of  fitting  a nonlinear  curve  to  a set  of  points  consists 
of  the  following  steps: 

1.  Estimate  the  M parameters  { A } . 
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A- 2 


Calculate  t he  elements  of  it!1  ami  A,  n.in  I • j. . i i 


and  (b). 

3.  Calculatt  tlu>  column  mattix  of  improved  v 1 • • ■ 

parameters  , using  Equation  (11). 

4.  Calculate  using  Equation  (2). 

3.  If  ' exceeds  an  arbitrarily  established  convergent*  rit.-ri.n, 
treat  the  latest  value  of  *,A  » as  {A(  • and  repeat  Steps  thr.aj.  l.  ,. 

If  ' does  not  exceed  the  convergence  criterion,  then  tin  lat.  ti  t. 
the  M parameters  ■'  \ } are  the  ones  desired  f r a best  tit,  mi  t'  < ,t 

squares  sense. 
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Mathematical  Description  of  Fiver  Plate  Motion 


The  posl t lon-t ime  response  of  a flyer  plate  element  is  characterized 
by  a smooth  curve  with  zero  initial  displacement  and  velocity,  approaching 
terminal  velocity  at  late  time.  The  rate  at  which  the  position-time  curve- 
approaches  the  terminal  velocity  asymptote  varies  from  shot  to  shot. 

Inspection  of  the  graph  of  a hyperbola  in  any  standard  text  on 
analytic  geometry  shows  that  an  ordinary  (second  order)  hyperbola  satisfies 
the  conditions  of  zero  initial  displacement  (with  an  upward  shift  of  the 
origin),  zero  initial  slope  and  terminal  asymptote.  All  that  is  needed 
in  addition  is  to  introduce  a parameter  which  will  control  the  rate  at 
which  the  curve  approaches  the  terminal  asymptote.  This  is  accomplished 
by  letting  the  exponent  be  a variable  parameter,  so  that  the  basic  form 
of  the  equation  tor  an  nt*1  order  hyperbola  becomes 


To  satisfy  the  condition  of  zero  initial  displacement,  set 


y = - b 
o 


(2) 


and  for  mathematical  convenience  in  the  derivations  to  follow,  assume 
that 


t = 0 
o 


and  set 


(3) 
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F.quation  (B-l)  then  takes  the  form 
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or,  solving  for  y. 


y ■ b i 1 + (at) 
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Figure  H-I  illustrates  how  the  parameters  of  an  n ' order  hyp«rt>>’ 
can  he  determined  graph  lea  1 1 v using  the  terminal  asytnpt  >te.  lie  • minil 
asvmptote  intersects  the  r axis  at  t ' h,  and  inters,  ■■  ('•  V i r 

v = -!>.  The  ordinate  to  the  curve  at  the  point  where  the  termini!  asv"pt  « 
intersects  the  T axis  is  y = V (2*  - l)h,  so  that 

n = JPJLl (7) 

iog(b  + 0 


The  velocity  at  any  time  is  given  by  the  equation 

n-  1 
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and  the  acceleration  is  given  hv  the  equation 


y’  = a~b  ( n—  1 ) 
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where 


p * instantaneous  detonation  pressure 
p * flyer  plate  mass  density 
H * flyer  plate  thickness 
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Inspection  of  Equations  (6),  (H),  and  (9)  shows  that 
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The  time  at  which  t tie  flyer  plate  displacement  equals  the  standoff 
distance,  s,  is  the  time  of  impact,  t.,  where 
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The  impact  velocity  is  therefore 


and  the  collision  angle  is 


(lh) 


(17) 


Note  that  this  entire  analysis  assumes  the  motion  of  the  flyer  plate 
to  he  purely  vertical.  When  n -2  the  peak  acceleration  occurs  at 


1 


.mil  the  peak  acceleration  is  given  by  the  expression 
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Because  t he  theoretical  acceleration  ! • infinite  u t ) Cot  a . . , a it i 

1.18)  and  (19)  cannot  be  used  to  calculate  rise  time  and  peak  -let*  at 
t ion  pressure  for  all  conditions.  Had  i double  exponential  ex(>l<  -!vi 
pressure  pulse  been- assumed,  e.g.  | I AO  IBS  KN  AND  AYKE  (1938,  132,  1 1 i 
1-11)]  this  particular  problem  would  have  been  avoided.  Mils  tvpi  ■ 
curve  is  employed  at  the1  Air  Force  Weapons  I.aboratorv  to  describe  i 
electromagnetic  pulse  (I'.MP)  time  history. 

In  obtaining  the  parametric  derivatives  needed  tor  least  |uur>  • 
curve  fitting,  use  was  made  of  the  fact  that  if 

Y « UV  (20) 

which  can  be  written  in  the  form 
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The  equations  for  the  required  partial  derivatives  are  as  follows: 
In (at)  In 
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Sine  i*  v is  a continuous  i unction,  cross  derivatives  are  o<hm1. 
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The  following  section  constitutes  i complete  description  >f  tin 
electronic  inst  rument.it  ton  .yster  used  it  the  USAI A l.xploive 
Impulse  holding  Research  Site. 

A.  The  I IH  Research  conJuctcd  at  this  facility  required  th 
measurement  of  various  data  including.  !>uf  not  mc<  irilv  liratt  • 
to,  explosive  detonation  pressure,  interface  plate  collision 
pressure,  detonation  velocity,  and  a flash  radiograph  of  the 
colliding  plates.  Several  unique  problems  had  to  h.  n-sol  >.a- 

1.  All  the  data  for  each  test  must  be  acquired  during 
the  time  >f  burning  ol  a slow  detonation  velocity  explosive 
with  detonation  velocity  in  the  range  Vj  - 2S00  m/s  - 4000  m/s. 

Th  i s represents  a total  time1  of  between  about  'JO  - 1 38  , sec , assum- 
ing a 12"  long  specimen  and  10  .sec  burning  time  for  the  line 
wave  generator. 

2.  Triggering  circuits  must  be  designed  so  that  all  of 
I lie  measurement  equipment  will  he  reliably  triggered  at  the 
instant  the  blasting  cap  initiates  the  detonation. 

3.  The  proper  time  relationships  between  the  various 
measuring  devices  must  be  accurately  maintained. 

4.  The  data  obtained  must  be  permanently  and  accurately 
recorded  for  later  analysis. 

5.  The  radiographic  film  cassette  and  mounting  terminal 
strip  for  the  wiring  to  the  specimen  must  he  protected  from  the 
effects  of  the  blast  and  shock  generated. 

C-l 
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B.  Sol ut ions  to  the  Problems 


1.  The  triggering  circuits  were  a two-port  problem. 

The  first  was  a method  suitable  for  utilizing  n "make"  type  f 
switch  to  trigger  the  electronic  at  time  I 0.  The  second 

was  generating  a pulse  of  the  proper  amplitude  and  length  to  cause 
reliable  triggering  of  the  equipment.  Ihe  flash  X-Ray  system 
used  in  conjunction  with  the  stress  measurement  sy.--.tcm  required  a 
triggering  pulse  of  a certain  amplitude  and  duration  .md  this 
pulse  was  also  used  to  trigger  the  other  electronics. 

2.  A "make"  type  of  switch  was  constructed  from  . f 1 1 > • 
thick  copper  foil  and  thin  pasteboard,  and  placed  directly  beneath 
the  blasting  cap.  Details  of  this  switch  are  shown  in  figure  C-l. 
When  the  cap  fired,  the  upper  foil  was  forced  through  the  hole 

111  the  pasteboard  and  contacted  the  lower  foil,  thereby  completing 
the  circuit. 

3.  A "trigger"  bo*  was  constructed  to  provide  the  necessar) 
pulse  output,  and  was  connected  to  the  "makd' switch , (herealter 
referred  to  as  the  "blast"  switch)  b>  means  of  #18  AW'D  2 conductor 
shielded  wire.  Details  of  the  trigger  box  are  shown  in  figure  ■.  - 

A pushbutton  switch  on  this  box,  essentially  in  parallel  with  t he- 
ll last  switch  was  used  to  trigger  the  electronics  for  set-up, 
prior  to  the  actual  shot. 

4.  The  output  of  the  trigger  box  was  transmitted  through 
KG58/U  coaxial  cable,  connecting  the  box  to  the  inputs  of 
Oscilloscope  #1,  (Time  base  #1  trigger  input).  Oscilloscope  #3 
(Time  base  #1  trigger  input),  Flash  x-RayTime  Delay  Input  Amplifier, 
Channel  #1  of  the  Pulsar  Power  Supply,  and  Channel  #1  of  the 
Abtronics  Time  Delay  Pulse  Generator.  (Channel  #2  of  this  unit  is 
not  used).  Additionally,  the  output  of  channel  #1  of  the  Abtronics 
was  connected  to  the  input  of  Channel  #2  of  the  Power  Supply  and 
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to  the  time  base  ”1  trigger  input  of  Oscilloscope  s 2 . >ii><  • • p< 

■ J anil  (Channel  !>2  of  the  Power  Supply  were  used  for  the  interta  < 
pressure  measurement,  the  appropriate  time  delay  ft  >n  t i mi  1 
could  he  dialed  irto  the  Ahtronics.  It  wa-  necessary  to  d'lav  ' 

Power  Supply  output  for  this  measurement  because  the  pul  ”;t  ;•  ■* 
has  a maximum  length  of  about  100  usee,  and  the  interface  press  ir* 
pulse  could  arrive  at  the  gage  anywhere  from  about  SO  usee  to  <■>.  i 
100  .sec,  depending  on  the  detonation  velocity. 

5.  Oscilloscope  8 I was  used  for  the  detonation  pri  sure 

measurement.  The  power  supply  pulse  was  observed  on  the  bottom 
trace  (using  time  base  "A")  and  the  magnified  pressure  pulse 
(superimposed  on  top  of  the  power  supply  pulse)  was  observed  on  the 
top  trace  (using  time  base  "B") . A type  "D"  "L";  "A”,  "0A", 

"53/54" , or  similar  iype  of  preamplifier  plug-in  could  be  used  for 
the  bottom  trace,  and  a type  "W"  or  "2"  could  be  used  for  the  top 
trace.  The  "K"  typo  ol'  preamp  was  found  to  he  entirely  Sat  i sfac*  or 

6.  Oscilloscope  "2  was  used  for  the  interface  pressure 
measurement.  The  configuration  here  was  the  same  as  for  Oscillo  ijw 
8 1 . except  that  the  input  triggering  pulse  was  delayed  by  tin 
amount  of  time  dialed  into  the  time-delay  generator  (from  time 

I = 0) . Additionally,  the  top  trace  of  Scope  « I always  needed 
some  time  delay  (with  respect  to  the  bottom  trace,  or  1 - 0) , 
because  of  the  difference  in  the  horizontal  magnification  factor. 

It  also  sometimes  proved  convenient  to  add  some  delay  to  the  top 
trace  of  Scope  *2  depending  on  several  factors.  These  factors, 
as  well  as  the  complete  operation  of  the  Oscilloscope  and  circuits 
will  be  described  later. 

7.  Oscilloscope  *3  was  used  for  measurement  of  the  Detonation 
Velocity.  It  was  connected  to  a specially  constructed  Detonation 
Velocity  probe  (glued  to  the  specimen)  by  means  of  RG5K/U  cable 

and  to  a specially  constructed  constant  current  source  that  provided 
a constant  current  output  for  the  probe.  Details  of  this  circuit 

04 


configuration  .ire  shown  in  figure  C-2.  lor  additional  information 
on  the  Flash  X-Ray  System  concerning  troubleshooting,  repair 
calibration,  etc.,  consult  the  Instruction  Manual  Model  730/2650 
(Fcxitron  Pulsed  X-Ray  Systems),  readily  available  from  Field 
(mission  Corp. , Melrose  Ave.  at  Linke  St.,  McMinnville,  Oregon, 

9"128. 

('.  Setup,  Calibration,  and  Operation  of  the  Complete  System 
for  a Fully  Instrumented  Test. 

1.  Open  both  bottles  of  gas  and  adjust  system  Freon 
chamber  pressure  for  about  19  PSI,  system  tubehead  pressure  for 
18  - 20  ounces;  and  nitrogen  pressure  for  25  PSI. 

2.  Activate  X-Ray  Console  main  power  switch,  and  make 
certain  High  Voltage  key  switch  is  in  "OFF"  position.  Ascertain 
that  "power  on"  lamps  are  lighted  for  Abtronics  time  delay 
generator.  Pulsar  Power  On  supply,  and  X-Ray  Delay  Trigger  amplifier. 
For  above  items  #1  and  #2  refer  to  Figures  C-5a  and  C-5b  for  item 
locat ion . 

3.  Turn  on  the  three  oscilloscopes,  allow  to  warm  up  for 
at  least  20  minutes  and  make  sure  constant  current  box  is  turned 
"OFF".  Insure  that  all  BMC  cables  are  connected  as  per  Figure  C-4. 

4.  Prepare  specimen  to  be  welded  as  indicated  in  Figure  C-6. 
Ascertain  velocity  probe  measures  - 95  - 105  ft  before  and  after 
being  glued  down  on  specimen  and  after  copper  ground  strip  is 

taped  to  probe.  Take  specimen  to  detonation  arena  (after  being 
taped  together),  place  in  position  for  test,  and  measure  distance 
from  end  of  implanted  gages  to  terminal  strip.  Take  specimen 
back,  cut  wires  tnis  distance,  strip  insulation,  and  solder  to 
appropriate  terminal  ends  on  gages.  Sec  Figures  C-4  and  C-6  for 
details  and  proper  polarity. 
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Blast  Switch  Trigger  Box  All  Coax  = RG58/U 
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FIGURE  C-6a 

ENTIRE  ASSEMBLY  TAP  HI)  TOC  KTTI  KR 
(not  shown) 


Pressure  Gage 
Leads  (Scope  #1) 


FIGURE  C-6  WIRING  DETAILS 
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:».  Place  prepare*!  speciRiei  in  proper  location  for  t»-  t. 
Ihts  ;» i 1 1 normally  lie  on  opposir  ->  ldi  of  steel  anvil  *i  ■ V Pay 
tuhehead,  and  somewhat  oil  n r ■ r t the  nrht  i • i • .*  tl  i id 

to  allow  the  le  laved  X-Ray  pulse  •*.  *.  • • . t r upon  t!  • a.  I at< 

All  tests  must  be  placed  upon  one  <•*  the  expendable  an  i!  • 
protect  the  largo  anvil.  On. a n • I ;u  non  . < nt.i  ly,  t« . r.  ! ■ 
remove  the  tuhehead  cover  p.  ite  and  look  between  the  ; lar  at 
the  1*1  on  the  tuhehead.  Hi  is  must  be  plain!)  o ihli  r i 
proper  X-Kiy,  and  if  it  is  not,  adjust  the  tuhehead  . t h. 
down  accordingly,  using  piece-  o:  styrofoam.  RI.I’LACI.  1 ' 'Hl.lir V • M 

<>.  Connect  gage  wires  t • terminal  strip,  observing  prop- 
polarity.  See  figure  i-7  tor  details.  Connect  blast  -witch  ■ 
per  figure  C -7,  and  short  terminals  together  a tew  tins  to  mil., 
certain  triggering  circuits  are  l).k.  Other  person  turn  triggei 
box  "OV  and  observes  yellow  "High  Voltage"  lamp  on  liel  .ved 
trigger  Xmplifier.  (Sec  Figure  C-5).  Lamp  will  Mink  ev-ry  time 
terminals  on  switch  are  shorted  (If  there  are  no  problems). 

Connect  Velocity  Probe  to  .ingle  COAX  line.  Miter  « ■ re 
to  center  conductor,  ground  to  ground  as  per  figure  i-7.  \ • 

alligator  clip  is  soldered  to  a tew'  inches  of  #28  wire  which 
then  soldered  to  the  coax  center  conductor.  This  clip  i then 
carefully  attached  to  the  probe  center  wire,  and  the  assembly 
taped  in  place.  Rl  SUR1:  no  pressure  is  exerted  on  the  probe,  or 
shorts  nr  "opens"  will  occur. 

8.  Using  the  black  cloth  film  "change  bag",  load  the  film 
holder  with  one  sheet  of  X-Ray  film.  Then  place  the  holder  in 
the  armored  film  cassette  using  the  s/8"  nuts  anJ  bolt  provided, 
making  certain  the  (»)  on  the  film  holder  faces  the  front  o!  the 

cassette  (the  1/8"  expendable  aluminum  face  piece).  Place  the 
cassette  in  position  just  behind  the  specimen  to  center  tl-  X-Uay 
photo  as  much  as  possible. 
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9.  Now  load  the  explosive  with  line  wave  generator  attach- 
ed. (See  Figure  C-8).  Tape  down  the  blast  switch  on  a wooden 
block,  and  tape  the  apex  of  the.  line  wave  generator  directly 
over  it.  This  is  where  the  blasting  cap  will  be  located. 

10.  At  this  point,  set  up  of  the  experiment  is  complete; 
setup  and  calibration  of  the  electronics  now  begins.  All 
scope  and  X-Ray  settings  should  be  entered  on  the  check-off  sheet 
provided,  for  complete  data  analysis  after  the  test  (Figure  C-ll). 
Refer  to  Figure'  C-9: 

a.  Set  Scope  «1  Time  Base  21A  to  single  sweep  mode. 

b.  Set  Time  Base  21A  for  lower  trace  and  Time  Base 
22A  for  upper  trace.  Set  Time  Base  22A  to  "SWEEPS  ONCE  FOR  EACH 
A DEL'D  B TRIG". 

c.  Set  both  "W"  and  "L"  plug  in  amps  to  "DC  Coupled" 
(Input  Attenuation  Control). 

d.  Set  "W"  controls  as  follows: 

Display  = A - V Comparison  Voltage  = 0.000 

V RANGE  = +11  Millivolts/cm  = 50 

c 

e.  Set  "L"  controls  as  follows: 

INPUT  = "AC" 

VOLTS/cm  = 1 V/cm 

f.  Set  Time  Base  #21A  controls  as  follows: 

SLOPE  = + SWEEP  = SINGLE  SWEEP 

COUPLING  = AC  TIME/cm  = 10  p volts/cm 

SOURCE  = EXT 


t 
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Explosive  set  on  top  of  specimen,  then 
blast  switch  taped  to  wood  block  and 


ERIC  L.I.VI.I.  - set  sn  reset 

lamp  remains  lit  until  deliberately  triggered 

g.  Set  I line  liase  :'22.\  controls  as  follows: 

SLOPE  - * 

COUPLING  •=  At! 

SOIJRCI  LOWER  PLUG  IN 
i 1 Ml  /cm  = 5 I-  volts/ cn. 

SKI. I P SKI  EPS  ONCt  FOR  EACH  A bid. 'VI)  B TRIG 
TRIG  LEVI  I.  Same  as  for  #21 

11.  Turn  trigger  box  "ON".  Every  time  push  button  is 
depressed  both  top  and  bottom  scope  traces  will  sweep.  (After 
each  sweep  reset  button  on  scope  must  be  depressed  and  will  stay 
lighted  until  trace  is  triggered) 

12.  temporarily  change  "input  Atten.  "to  "GND".  Now 
operate  the  "Intensity",  "Focus",  "Astigmatism",  and  "Vertical" 
and  "Horizontal"  position  controls  to  place  the  bottom  trace 
(dot)  on  the  bottom  line  of  the  graticule  at  the  left-hand  edge 
("L"  plug  in).  Now  operate  similar  controls  relative  to  the 
"K"  plug  in,  and  place  the  trace  (dot)  on  the  4th  line  down  from 
the  top  and  at  the  left-hand  edge.  Check  these  adjustments  several 
times  to  be  sure  there  is  no  appreciable  drift  especially  with  the 
top  trace.  Screen  illumination  and  top  trace  intensity  are  very 
important  and  must  be  carefully  balanced  so  that  the  pressure 

peak  will  have  as  much  contrast  as  possible.  Pencil  marks  have 
been  made  next  to  the  Screen  illumination  and  trace  intensity 
controls  for  optimum  settings  and  should  lie  followed. 

13.  Now,  change  "Input  Atten."  to  "DC"  and  operate  the 
scope  reset  button;  then  operate  the  trigger  box  push  button  while 
observing  the  bottom  trace.  A "step"  similar  to  that  of  Figure 
#9a  should  be  observed.  At  this  point,  using  a small  screwdriver, 
adjust  the  pulse  width  control  on  channel  #1  of  the  Pulsar  Power 
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Supply  (Sec  Figure  C-5)  to  obtain  a pulse  length  of  about  7 cm 
(70  usecl.  Also,  Jo  not  pulse  the  Power  Supply  more  than  about 
six  times  per  minute. 

II.  Move  the  "Input  At  ten"  setting  on  the  "W"  plug-in 
back  to  the  "DC"  position.  After  observing  about  how  high  the  "step 
on  the  bottom  trace  is  (usually  about  1.200V),  dial  that  much  in 
on  the  Comparison  Voltage  dial  (V^J.  This  should  bring  the  top 
of  the  step  onto  the  scope  face  (TOP  TRACT).  Then  carefully  auimt 
for  more  or  less  V to  bring  the  top  trace  down  as  close  as  pos> 
lble  to  the  top  trace  base  line  (4th  line  down  from  top).  Hnter 
these  values  as  noted  on  the  test  checkoff  sheet. 

15.  Set  the  time  delay  vernier  dial  for  the  proper  amount 
of  time  delay.  Tacit  major  division  corresponds  to  10  microseconds 
at  this  sweep  rate  (10  psec/cm  on  the  BOTTOM  TRACI. ) . Since  there 
is  an  effective  50  usee  "window"  on  the  top  trace,  the  time  delay- 
should  be  adjusted  to  try  to  center  the  expected  detonation  stress 
pulse  in  this  "window".  To  the  total  amount  of  time  for  the 
explosive  to  burn  to  the  top  gage  must  be  added  about  16  psec/inch, 
burning  time  for  the  line  wave  generator.  For  an  average  value  of 
5000  m/ s detonation  velocity  one  can  calculate  about  10  sec/ inch, 
to  which  must  always  be  added  the  16  usee  line  wave  generator 
burning  time.  If  the  detonation  pressure  gage  is  located  4 in. 
down  the  plate,  for  example,  one  could  expect  the  pulse  to  arrive 

it  about  56  usee  from  time  T = 0.  To  center  the  top  trace  "window" 
at  about  55  or  56  usee , one  would  need  to  dial  in  about  50  usee 
time  delay  into  scope  #1.  For  a typical  result  see  Figure  C-lOc. 

16.  Insure  that  both  inputs  to  Scope  #2  are  in  "CND" 
position,  ami  proceed  to  "baseline"  noth  traces  in  tne  same  maimer 
as  with  Scope  #1.  Once  baselined  and  stable,  set  inputs  to  "IK 
position  and  locate  top  trace  in  the  same  manner  as  with  Scope  «1. 
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TOP  TRACE 


FIGURE  C-10 
OSCILLOSCOPES  #1  ft  #2 


VERT  = 50  mV/cm 
HORIZ  = 5 usec/cm 

BOTTOM  TRACI: 

VERT  = 1 V/cm 
HORIZ  = 10  usec/cm 
PULSE 

WIDTH  - 70.2  usee 
PULSE 

HEIGHT  e 1.2  V 
TIME 

DELAY  = 30  usec 
(TOP  TRACE) 


HGIIRI  C - ] ()a  - SET  UP  FOR  BOTTOM  TRACE 


TOP  TRACE 

VERT  = 50  mV/cm 
HORIZ  = 5 usec/cm 
BASELINE 
DIFF  = 0.3  cm 

BOTTOM  TRACE 

VERT  = lV/cm 
HORIZ  = 10  usec/cm 
PULSE 

WIDTH  = 70.2  usec 
PULSE 

HEIGHT  = 1.2V 
TIME 

DELAY  = 30  usec 
(TOP  TRACE) 

VOLTAGE  OFFSET 
(Vc)  = 1.185  v 


HGURI  C- 1 01)  - SI  I UP  FOR  TOP  TRACI 


C - 2 5 


FIGURE  C-10 

OSCILLOSCOPES  *1  6 *2  (Contd) 


TOP  TRACE 


VERT  =50  mV/ cm 
HORIZ  = 5 ysec/cm 

BOTTOM  TRACE 

VERT  = 1 V/cm 
HORIZ  = 10  ysec/cm 
PULSE 

WIDTH  = 80  ysec 
PULSE 

HEIGHT  = 1.25  V 

V = 1.235  v 

c 

TIME 

DELAY  = 40  ysec 
(TOP  TRACE) 


FIGURE  C - 10c  - DETONATION  PRESSURE  PULSE  PHOTO 
(PHOTO  WAS  P RE-BASELINED) 


TOP  TRACE 

VERT  = 50  mV/cm 
HORIZ  = 5 ysec/cm 

BOTTOM  TRACE 

VERT  = 1 V/cm 
HORIZ  = 10  ysec/cm 
PULSE 

WIDTH  = 80.4  ysec 
PULSE 

HEIGHT  = 1.4V 

V = 1.384  V 
c 

TIME 

DELAY  = 2 ysec 
(TOP  TRACE) 

TIME  DELAY  (FROM  T = 0 
BY  TIME  DELAY  GEN  = 80  ysec 


FIGURI  i - Kiel  - INTERFACE  PRESSURE  PULSE  PHOTO 
(PHOTO  WAS  PRE -BASELINED) 
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FIGURI  C-l  I 


FINAL  CHI  l!KL  1 SI  FOR  QUADRUPLI 
INSTRUMENT  SHOT 

USAFA  EXPERIMENT  * 


A.  Before  assembling  specimen  on  firing  table,  in  ir  that 
each  trigger  line  triggers  its  respective  device  at  hast 
10  times  consecutively. 

B.  Inspect  all  cables  for  cuts  and  bruises,  and  insure 
maximum  separation  from  Flash  X-Ray  leads.  Insure 
proper  cable  connections  at  arena. 

C.  After  assembling  specimen  on  firing  table,  insure 
that  specimen  is  properly  positioned  for  X-Ray  shot. 

Insure  that  "X"  on  tube  cover  face  cun  be  seen  by  looking 
between  plates  of  specimen. 

D.  Install  assembled  armored  cassette  in  proper  position. 

F.  INSTALL  PROTLCTIVL  COVER  PLATE  OVER  TUBE  HEAD! 

F.  Have  initiation  device  connected;  FIRING  CIRCUl'l 
"READY”;  and  firing  switch  close  at  hand  in  instrumentation 
faci 1 ity . 


CHECK  LIST 

1.  Insure  three  scope  cameras  arc  operational  and  contain 
f i lm. 

2.  Insure  flash  X-Ray  gas  pressures  are  correct  and  that 
trigger  lamp  and  hi -voltage  lamp  are  lit. 

3.  Inspect  all  settings  on  scope  #1  for  correctness. 

Insure  time  delay  is  proper  value.  Insure  proper 
connections . 

4.  Pulse  first  channel  of  Pulsor  with  trigger  box.  insure 
that  pulse  appears  on  bottom  trace.  Zero  both  traces 
several  times  to  insure  minimum  baseline  error. 

5.  Repeat  triggering  of  Pulsor  and  dial  in  offset 
voltage  on  "W"  plug-in  amp  until  top  of  pulse  appears  in 
top  window.  Carefully  adjust  until  top  of  pulse  is  as 
close  as  practicable  to  upper  trace  baseline.  At  this  time, 
enter  the  following  values: 
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FIGURE  C 11  (Contd ) 


UPPER  I'RACi: 


A.  VERTICAL  SCALE: 

R.  HORIZONTAL  SCALE: 

C.  OEESET  VOLTAGE  AS  INDICATED  ON  DIAL: 

I).  DISTANCE  BETWEEN  TOP  01  PULSE  AND  UPPER  I' RACE 
BASELINE : 

E . IT  MI  DELAY  AS  INDICATED  ON  DIAL: 

LOWER  TRACE 


G.  VERTICAL  SCALE: 

H.  HORIZONTAL  SCALE: 

I.  HEIGHT  OF  PULSE: 

LENGTH  OF'  PULSE: 

K.  TIME  DELAY  FROM  TIME  ZERO  IF  USED: 

L.  DISTANCE  OF  PRESSURE  GUAGE  FROM  INITIATION  INI)  0‘ 
SPECIMEN: 

0.  Inspect  all  settings  on  Scope  #2  for  correctness. 

Insure  time  delay  is  proper  value.  Insure  proper  connect i ms. 

7.  Pulse  second  channel  of  Pulsor  with  trigger  box.  insure 
that  pulse  appears  on  bottom  trace.  Zero  both  traces 
several  times  to  insure  minimum  baseline  error. 

8.  Repeat  triggering  of  Pulsor  and  dial  in  offset 
voltage  on  "W"  pin  -in  amp  until  top  of  pulse  appears 
in  top  window.  (aref.Ily  adjust  until  top  of  pulse  is 

as  close  as  practicable  to  upper  trace  baseline.  At  this 
time,  enter  the  fol  Owing  values: 


upper  trace: 


A.  VERTICAL  SCALE: 

B.  HORIZONTAL  SCALE: 

C.  OFFSET  VOLTAGE  AS  INDICATED  ON  DIAL: 
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ITGIIRI  (Ml  (Conti!) 


I).  DISTANC'  BETWEEN  IOP  OF  Ptll.SI  AND  UPPER  TKA(.| 
BASl  LINE: 

I..  MM!  DELAY  AS  INDICATED  ON  DIM.: 


low i k nun. 


F.  VLRT1CAL  SCAM  : 

G.  HORIZONTAL  SCALL: 

H.  HEIGHT  OF  PULSE: 

I.  LENGTH  OF  PULSI  : 

.).  DISTANC!  OF  INTERFACE  PRESSURE  GUAG!  FROM  INITI 
AT  ION  END  OF  SPECIMEN: 

NOTE.:  ABTRONICS  TIME  DELAY 

9.  Inspect  all  settings  on  Scope  *3  for  correctness.  Insure 
proper  connections. 

A.  VELOCITY  PROBE  CABLE  TO  "PROBE"  ON  CONSTANT 
CURRENT  BOX 

B.  "SCOPE"  OUTPUT  TO  SCOPE  INPUT 

c.  scope:  set  on  external  trigger 

D.  SET  VERT . SCALE  TO  2V/(  AND  HORIZONTAL  SCALL  10 
10  M SEC/CM.  11  ANOTHER  VALUE  FOR  HORIZONTAL 
SETTING,  ENTER  HERE 


10.  With  Box  Switch  in  "OFF"  position,  baseline  scope 
trace  several  times  to  insure  minimum  baseline  error. 

11.  MEASURE  VELOCITY  PROBE  CABLE  BETWEEN  CENTER  AND 
GROUND.  SHOULD  BE  IN  VICINITY  01  100- ISO  OHMS.  I E 
MARKEDLY  DIFFERENT  FROM  THIS,  RESIT  FIRING  CIRCUIT  SAFETY 
SWITCHES,  AND  REINSPECT  PROBE  AT  SPECIMEN,  AND  REPAIR 

AS  NECESSARY. 

12.  Turn  box  switch  to  "ON",  setting  to  "CA1 IRRATE";  and 
adjust  potent io-meter  so  that  the  trace  moves  up  exactly 
4 CM  from  baseline  when  the  switch  is  turned  on.  INSURE 
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I- 1 <;UKI  C-ll  (Conti!) 


I'HAT  RESII  LAMP  REMAINS  LIT  IIVI1I.  TRIGGERED  lit  TRIGGIK 
BOX . 

IT.  Turn  switch  "OFT";  turn  setting  to  "PROBE";  react iva* 
switch  and  insure  that  trace  moves  up  at  approximate!)  tin 
1 CM,  depending  on  the  length  of  the  prohe. 

14.  Return  to  X-Ray  console;  recheck  both  gas  pre-  ires 
for  correctness;  time  delay  for  proper  setting;  and  insure 
again  that  trigger  and  High  Voltage  lamps  arc  lit.  In?  • 
the  following  at  this  time: 

A.  FRLON  I’RLSSURI  : 

B.  NITROGEN  PRFSSUR! 
t:.  TIM!  DELAY  SETTING: 

15.  Return  to  scope  1 ; recheck  both  traces  for  proper 
operation.  Slake  sure  bottom  trace  brightness  is  turned 
down  just  out  of  sight.  Insure  that  screen  illumination  is 
turned  down  to  the  mark  on  the  scope.  Insure  "RESET" 
button  is  lit. 

lb.  CLOSE  0 1 CAMERA  VIEWING  PORI! 

17.  Return  to  Scope  #2;  recheck  both  traces  for  proper 
operation.  Make  sure  bottom  trace  brightness  is  turned 
down  just  out  of  sight.  Insure  that  screen  illumination 
is  turned  down  to  the  mark  on  the  scope.  Insure  "RESIT" 
but  ton  is  lit. 

18.  CLOSE:  * 2 CAMERA  VIEWING  PORT! 

19.  Return  to  Scope  #3;  recheck  to  see  that  scope  triggers 
properly.  Make  sure  trace  brightness  is  turned  down  just 
out  of  sight.  Insure  that  screen  illumination  is  turned 
down  to  the  mark  on  the  scope. 

20.  CLOSE  »3  CAMERA  VIEWING  PORT! 

21.  BASELINE  SCOPES  II 1 & 112'. 

22.  Activate  High  Voltage  and  charge  X-Ray  to  proper 
voltage,  being  careful  to  charge  slowly  to  preclude 
prefi  res . 

23.  TURN  CONSTANT-CURRENT  BOX  SWITCH  "ON". 
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FIGURE  C-ll  (Contd) 


24.  ACT  I VAT!  "RESET"  BUTTON . 

25.  OPEN  #3  CAMERA  SHUTTER . 

26.  OPEN  "2  CAMERA  SHUTTER. 

27.  OPEN  «1  CAMERA  SHUTTER. 

28.  Glance  at  all  equipment ; insure  that  all  three  scope 
reset  buttons  are  lit,  and  that  X-Ray  is  charged  and  both 
lamps  on  trigger  amplifier  are  lit.  Insure  three  lamps 
on  pulsor  are  lit. 

29.  GIVE  WARNING  AN!)  PRESS  FIRING  CIRCUIT  SWITCH. 

30.  CLOSE  ALL  3 CAMERA  SHUTTERS . 

31.  Operate  red  button  on  #1  camera;  lift  cutter  bar 
slightly;  grasp  paper  firmly,  and  null  out  to  stop  -- 
discard  paper. 

32.  Repeat  step  »30  for  # 2 camera . 

33.  Repeat  step  ^30  for  0.3  camera. 

34.  Turn  high-voltage  dial  on  X-Ray  console  down  to  21. RO. 
Turn  "READY"  switch  on  trigger  amplifier  "OFF". 

35.  After  count  of  about  "15"  from  start  of  film  developing 
procedures,  open  cameras,  retrieve  photos,  and  coat  with 
fixing  chemical  provided. 

36.  When  photos  are  dry,  enter  appropriate  data  on  backs. 

37.  Perform  POST-SHOT  cleanup  tasks. 

.38.  TURN  OFF  BOTH  CHANNELS  OF  TRIGGER  BOX. 

39.  TURN  SWITCH  ON  CONSTANT  CURRENT  BOX  "OFF". 

SPECIMEN  DATA 

40.  THICKNESS 

41.  STANDOFF  DISTANCE 

42.  LOADING 
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FIGURE  C-ll  (Contd) 


43.  TYPE  OF  EXPLOSIVE 

44 . REMARKS 


) 
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17.  The  time  delay  for  Scope  »2  is  somewhat  different. 

The  time  delay  dial  has  been  set  at  perhaps  one  or  two  micro- 
seconds. (.Just  enough  to  insure  proper  operation  of  the  delay 
circuits),  and  all  the  delay  necessary  has  been  dialed  into  the 
Abtronics  time  delay  generator.  With  the  top  and  bottom  traces 
coming  on  at  essentially  the  same  time,  the  time  delay  generator 
can  be  set  to  center  the  interface  pressure  pulse  on  the  top 
trace  of  Scope  #2.  1'hc  interface  pressure  pulse  can  be  expected 
to  arrive  at  the  interface  gage  anywhere  from  20  to  perhaps  SO 
usee  later  than  the  time  needed  just  for  the  additional  distance. 

This  is  caused  by  the  plate  stand-off  distance,  plate  thickness, 
etc.  With  a distance  of  6 in  down  the  plate  for  the  interface 
gage,  using  1/4"  thick  6U61-T6  aluminum  plates,  the  interface  pulse 
has  been  appearing  at  about  80  to  00  usee  from  T = 0.  Therefore 

a good  delay  is  60  usee  for  the  time  delay  generator.  (Assuming 
3000  m/s  detonation  velocity). 

18.  Since  "spikes"  are  produced  on  the  power  supply  "step" 
pulse  when  the  time  dealy  generator,  power  supply,  and  Flash  X-Ray 
fire,  the  various  delay  times  need  to  additionally  reflect  some 
planning  in  keeping  the  areas  where  both  detonation  and  interface 
stress  pulses  are  expected  as  free  as  possible  from  any  interfer- 
ence. The  "spikes"  produced  from  these  instruments  are  of  a sufficient 
magnitude  to  obliterate  the  stress  pulses,  if  superimposed  thereon. 

The  Flash  X-Ray  produces  such  a large  interference  pattern  (A 

damped  oscillation  of  1 or  2 volts  peak-to-peak  and  perhaps 
lOpscc  long),  that  it  has  been  set  to  fire  after  both  stress  pulses 
have  been  recorded. 

19.  The  next  step  in  the  set-up  of  the  electronics  is 
Scope  * 3.  (Refer  to  Figure  C-12).  It  uses  only  the  top  trace 
(Plug  in  "L"  with  time  base  #21A)  and  is  triggered  at  T = 0 
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directly  from  the  trigger  box. 
TIMI:  BASH  21A  = UPPHR  BEAM 
SLOPE  = + 

COUPLING  = AC 
SOURCE  = EXT 
SWEEP  = SINGLE  SWEEP 


The  following  settings  are  made: 

T I MI:  /cm  - 10  usec/cm 
I RIGGER  LEVEL  = set  so  re  et  l imp 
remains  lit  until 
de 1 i berate  1 y 
triggered 

PLUG  IN  "D" 

INPUT  = DC 
VOLTS/cm  = 2 V/cm 


a.  With  the  constant  current  generator  still  in  the 
off  position,  and  using  the  Intensity,  Focus,  Astigmatism,  and  ver- 
tical and  horizontal  centering  controls,  center  the  top  trace  of 
Scope  #3  on  the  4th  line  down  from  the  top  of  the  graticule,  and  at 
the  left  hand  edge.  Now  put  the  constant -current  box  in  the  "CAL" 
position  and  turn  the  box  to  "ON".  The  trace  should  immediately 
rise  to  about  the  top  of  the  graticule.  Lsing  the  adjustment  control, 
set  the  trace  exactly  on  the  top  line  of  the  graticule.  Turn  the 
box  off;  the  trace  should  fall  back  to  the  4th  line  down.  Now  set 
the  box  to  "PROBE"  and  turn  on  again;  this  time  the  trace  will  rise 
to  whatever  value  of  resistance  the  probe  is.  Each  centimeter  of 
height  represents  2 volts  or  25  ohms.  The  resistance  wire  used  in 
the  probes  has  a value  of  86. 6 ohms/ft , so  a direct  relationship 
of  distance  to  time  is  established.  The  resistance  of  the  completed 
probe  plus  small  line  losses  is  usually  on  the  order  of  95  - 105 
ohms.  After  this  value  is  determined,  turn  the  constant-current 
box  off  and  enter  the  appropriate  readings  on  the  check-off  sheet. 

Be  sure  the  box  is  turned  OFF  until  just  before  the  event  because 
current  drain  on  the  battery  is  very  heavy  and  the  transistor 
will  also  get  very  hot  and  may  burn  up  if  left  on  for  any  extended 
length  of  time. 


C- 15 


wet  >-r  to 


20.  Now  cap  up  the  device,  and  ARM  the  tiring  circuit. 

Figure  C-13  for  the  cap  and  explosive  configuration. 

21.  Next,  dial  in  the  appropriate  amount  of  delay  on  the 
Flash  X-Ray  delayed  trigger  amplifier.  (See  Figure  C-5)  . Usually 
this  will  be  in  the  order  of  100  - 120  usee,  again  depending  on  the 
detonation  velocity.  This  must  always  be  set  to  fire  after  both 
stress  peaks  have  been  recorded  or  its  interference  may  destroy 
either  pulse.  lor  a detonation  velocity  of  SO 00  m/sec,  an  X-Ray 
setting  of  about  105  usee  would  probably  be  adequate.  Now  check  to 
see  that  gas  pressures  are  correct : 

FREON  = 19  lb 

NITROGEN  = 25  lb 

FREON  (TU BEHEAD)  18  - 20  ounces 

Make  sure  the  high-voltage  and  ready  lamps  are  lighted.  Turn  the 
high  voltage  key  switch  to  "ON",  and  depress  the  "HIGH  VOE'IAGE  ON" 
button.  The  red  lamp  should  come  on  and  one  can  now  charge  the 
system  SEOWLY  (Use  no  more  than  1/2  ma  charging  current)  to  whatever 
value  gives  good  X-Ray  resolution.  23  KV  appears  to  be  entirely 
adequate.  When  the  system  is  charged  check  the  following  items: 

a.  All  three  scope  reset  lamps  lighted. 

b.  All  three  scope  viewing  ports  shut. 

c.  Optionally,  before  the  Flash  X-Ray  System  is  charged 
up. 

Scopes  Ml  and  M 2 can  be  pre-base l i ned . This  is  accomplished  by 
opening  the  shutters  on  cameras  Ml  and  »2,  operating  the  trigger 
box  and  then  closing  shutters  Ml  and  M2.  This  provides  a baseline 
reference  against  which  the  actual  test  baseline  can  be  measured. 
This  measure  precludes  or  minimires  the  possibility  of  drift, 
especially  of  the  top  trace  (which  is  sensitive  to  temperature 
variations,  etc.). 
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22.  Connect  firing  switch  (w/battcry)  assy,  to  firing 
circuit  terminals  in  detonation  control  breaker.  (Run  through 
window  of  instrumentation  facility).  Check  again  visually  to 
make  sure  the  X-Ray  trigger  Amplifier  red  "RPAPY"  and  yellow  "IIY" 
lamps  are  both  lighted.  Insure  reset  buttons  on  all  three  . >pe 
are  lighted,  and  that  all  three  camera  viewports  are  closed. 

.13.  Open  all  three  camera  shutters  and  insure  constant 
current  box  is  in  "PROBF"  position  and  turn  box  "ON".  Give  warning 
on  speaker.  Now  press  firing  switch  to  detonate  device.  Turn 
constant  current  box  "OFF"  and  close  all  three  camera  shuLters. 

24.  Operate  the  red  button  on  the  back  of  each  earner..; 
grasp  end  of  paper  and  pull  straight  out  until  it  stops.  Now  count 
to  15,  then  open  the  door  on  the  back  of  the  camera  and  retrieve 
photo.  Coat  with  fi.xer  provided. 

25.  Turn  the  X-Ray  IIV  control  down  to  zero,  and  turn  the 
safety  key  switch  to  off. 

26.  Retrieve  the  specimen  and  armored  cassette.  f'he  film 
develops  very  nicely  in  one  of  the  90  - second  developing  machines 
for  medical  X-Ray  use. 

1).  Data  Analysis 

1.  For  both  the  detonation  and  interface  pressures,  the 
formula  to  find  the  pressure  in  kilobars  is 
V - V 

= •!— — 1 where:  P,  = Pressure,  in  Kilobars 

k V k 

Vp  = Height  of  the  pressure  pulse  above 
the  baseline  (MV) 

= Offset  baseline  (MV). 

V = Offset  voltage  as  indicated  on  the 
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Comparison  Voltage  Dial  (V  ) on  the  "W" 

c 

Plug  in  (VOLTS). 

If  the  pressure  in  P.S.I.  is  desired,  multiply  this  result  by  14.?. 

1.  For  the  detonation  velocity,  draw  a straight  line  to 
host  coincide  with  the  slope  of  the  detonation  velocity  photo. 

(See  Figure  C-14).  Since  the  vertical  scale  = 25  ohms 

per  cm,  and  the  known  resistance  of  the  probe  is  86.6  ohms/ ft  the 

velocity  is  equal  to  R /T,  where  R is  in  ohms  and  T is  in  usee  and 

V is  in  ft/psec.  0^(>c&urse,  to  express  the  result  Jn  m/ s , as  has 

^ ^ f)  1 

usually  been  the  case,  multiply  the  result  by  j 10  1 

W 

3.  A complete  record  of  all  Site  experiments  has  been  kept 
in  three  log  books.  One,  the  Master  Log,  contains  a record  of  all 
tests.  The  instrumented  shots  that  have  been  evaluated  with  a 
computer  are  located  in  another  log,  and  all  other  tests  are  in  the 
third  log  book.  Additionally,  a work  sheet  has  been  kept  for  each 
test.  Use  of  these  sheets,  especially  with  the  fully  instrumented 
tests,  tend  to  preclude  operator  error  due  to  the  relative  complexity 
of  the  system. 

L.  Additional  Information 

1.  The  velocity  probes  are  made  from  aluminum  tubing  and 
"skip"  wound  nylon  insulated  resistance  wire  with  a resistance  of 
8b. 6 ohms/ft.  The  tubing  has  an  inside  diameter  of  about  19  mils 
with  a wall  thickness  of  1 mil  and  is  accordingly  very  fragile. 

The  tube  is  measured  (usually,  one  ft  lengths  have  been  used) 
scored  with  a sharp  razor  blade  and  broken  off.  By  cutting  the  wire 
with  very  sharp  side  cutters  or  scissors  it  can  be  readily  threaded 
through  the  tubing.  One  end  is  then  pinched  shut  and  a drop  of 
"DlK.'O"  cement  or  equivalent  is  applied  to  the  other  end.  A small 
copper  strip,  with  two  or  three  inches  of  Insulated  wire 
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soldered  to  it,  is  placed  over  the  probe  at  the  connection  end  and 
taped  down  on  both  sides,  so  as  to  make  good  electrical  connection, 
liiir  not  w/rh  enough  pressure  to  col  lapse  the  tube  wall.  See 
I igure  t.'-l  - tor  completed  probe,  "'lore  information  on  this  subject 
is  contained  in  Document  OR- 70 1 1 -2  157  , (I  UV  DI.VLLOl’MbNT  PROGRAM 
PROGRESS  H1.P0R1,  Third  Progress  Review  Meeting,  dated  10/29/71. 

2.  Since  the  velocity  probe  is  so  fragile,  great  care 
must  be  taken  when  gluing  it  down  to  the  specimen  (using  the  same 
"IHX'O"  cement  or  equiv.J  and  when  attaching  the  coax  cable  when  the 
specimen  is  in  position  on  the  anvil.  Once  it  is  in  place,  remove 
the  COAX  fitting  at  the  constant -current  box  marked  "PROBP"  , and 
mea  ure  with  the  VOM.  The  reading  should  be  somewhere  in  t he  vicinity 
of  95  - 1()5.  If  markedly  different  from  these  values,  the  probe  i 
either  shorted  or  open,  and  appropriate  repair  is  necessary. 

.5.  Before  closing  the  site  for  the  day,  it  is  highly 
advisable  to  completely  shut  off  the  Freon  and  Nitrogen  bottles, 
because  both  systems  have  inherent  small  leaks  that  have  been  imposs- 
ible to  locate. 

4.  Before  running  a test  it  is  advisable  to  take  a test 
photograph  with  each  of  the  three  cameras  (with  the  specimen 
completely  wired),  to  insure  proper  operation  of  the  cameras  as 
well  as  all  the  systems  involved,  for  the  same  reasons,  it  is  advis- 
able to  "bleed"  bo*h  gas  systems  of  the  X-Ray  if  it  has  been  more 
than  a few  days  since  the  last  test.  A dosimeter  has  been  used  to 
measure  radiation  onput  from  the  X-Ray  at  various  intervals,  to 
insure  proper  operation.  Electric  blankets  were  placed  over 
both  pulsers  (located  in  a shed  next  to  the  detonation  arena) 
to  keep  the  freon  from  liquifying;  the  blanket  must  be  temporarily 
removed  from  the  pulser  in  operation  before  running  the  test 
or  the  high  voltage  pulse  discharge  will  find  a short-circuit 
across  the  blanket.  Additionally,  if  there  is  ANY  indication  that 
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FIGURE  C - 1 4 

oscilloscope  #3 


FIGURE  C-14a  - DETONATION  VELOCITY  TRACE 

(WITH  FULL  INSTRUMENTATION) 

..  R/86.6  75/  86 . 6 106 
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FIGURE  C - 1 4b  - DETONATION  VELOCITY  TRACE 
(Wini  FLASH  X-RAY  ONLY) 

v _ R/86.6  50/86.6  106 

0 T = 56  X 3.28  “ 


TOP  TRACE  ONLY 

VERT  = 2 V/cm 
HORIZ  = 10  ysec/cm 

NO  TIME  DELAY 


4,180  M/sec 

TOP  TRACE  ONLY 

VERT  = 2 V/cm 
HORIZ  = 10  u sec/cm 

Nf  DELAY 


3,140  M/sec 


C-41 


iaojM  ui 


the  Freon  has  liquefied  (such  as  attempting  to  bleed  the  Freon  and 
getting  any  result  other  than  a steady  hiss  or  the  body  of  the 
pulser  being  very  cold  to  the  touch),  do  NOT  pulse  the  system,  or 
serious  internal  damage  to  the  pulser  may  result. 

5.  With  reference  to  Figure  C-14,  the  extra  noise 
and  nonlinearity  that  seem  to  be  generated  in  13a  appear  to  be  a 
function  of  the  high  voltage  and  current  pulse  through  the  gages. 
Many  different  confi guartions  have  been  tried  in  an  attempt  to 
clear  up  this  problem,  but  to  no  avail.  Also,  during  velocity 
calculations  an  apparent  discrepancy  was  noted  between  the  detona- 
tion velocity,  as  calculated  from  the  oscillograph,  and  the  velocity 
as  calculated  from  known  distances  of  the  gages,  time  delay  of  the 
electronics,  and  time  that  the  pulses  arrived  at  the  gages.  This 
discrepancy  appeared  to  be  consistently  20  - 25  percent.  It  was 
finally  found  that  this  apparent  discrepancy  was  noted  during  the 
series  using  dynamite,  and  that  there  is  a large  difference  between 
the  initial  velocity  and  the  final  velocity  for  dynamite.  This 
difference  is  negligible  with  the  sheet  explosives,  as  is  evidenced 
by  many  oscillographs  taken  during  the  various  test  series.  In 
summary,  the  difference  between  the  initial  and  final  velocities 
ippeared  to  be  the  culprit,  because  average  velocity  calculations 
for  both  methods  closely  agree. 
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APPENDIX  D 


CHARACTERISTICS  OF  EXPI.OS  I VKS 

During  the  months  of  October  and  November  1969  it  was  decided 
to  repeat  some  explosion  welding  experiments  from  the  parallel  alu- 
minum welding  study  to  check  repeatability.  Forty  percent  extra 
dynamite  from  a different  lot  of  explosive  was  used.  Shear  and 
tensile  strength  tests  showed  the  weld  strengths  to  be  lower  than 
for  the  original  set  of  welds.  Additional  welding  experiments 
indicated  that  the  inferior  results  were  attributable  to  the  charac- 
teristics of  the  second  lot  of  explosive.  When  the  welding  xperi- 
ments  were  rerun  with  the  original  lot  of  explosive,  the  bond  strengths 
were  consistent  with  the  original  results.  It  was  eventually  deter- 
mined that  the  second  lot  of  explosive  contained  excessive  moisture 
associated  with  the  ammonium  nitrate,  an  ingredient  common  to  most 
medium  and  low  detonation  velocity  dynamites. 

These  results  cast  considerable  doubt  about  the  suitability  of 
granular,  nitrate-containing  dynamites  either  for  explosion  welding 
field  applications,  or  for  initiating  other  parametric  investi- 
gations. An  effort  was  begun  to  bring  into  the  program  another 
explosive  with  more  favorable  characteristics.  IRECO  Chemical 
Company  of  Salt  Lake  City,  Utah,  produces  a series  of  slurry  and 
gel  explosives  which  appeared  to  be  a reasonable  basis  for  developing 
an  explosion  welding  explosive.  After  some  discussion,  IRECO  officials 
agreed  to  modify  their  slurry  explosive  and  supply  enough  explosive 
to  complete  a series  of  preliminary  welding  experiments.  The  first 


D-l 


explosive  composition  supplied  proved  to  have  undesirable  physical 
characteristics:  viz,  it  was  a putty-like  composition, "sticky"  and 
difficult  to  apply  uniformly.  A second  explosive  composition  wa  ; 
supplied  by  IRF.CO,  and  on  the  strength  of  its  welding  performance  wa 
accepted,  although  the  physical  characteristics  were  not  considered 
to  be  ideal.  The  second  IRECO  explosive  was  a water  gel  of  proprietary 
composition.  It  had  a physical  consistency  like  "fresh  liver".  Its 
detonation  velocity  was  5000  m/sec  at  a density  of  1.5  g/cm\  and  it 
had  a relatively  high  impulse  value.  This  explosive  was  purchased 
in  sufficient  quantity  to  conduct  a steel-steel  preset  angle  weld- 
ing parametric  study. 

The  IRECO  gel  explosive  was  supplied  in  cast  12  inch  cubes.  In 
order  to  cut  sheets  from  these  blocks  as  thin  as  1/4  inch,  it  was 
necessary  to  develop  a special  cutting  fixture.  After  evaluating 
several  suggested  methods,  a "mitre  box  and  meat  saw"  technique  was 
adopted  with  successful  results.  A wooden  mitre  box  with  an  inter- 
nal sliding  tray  was  built.  The  mitre  box  supported  five  sides  ot 
the  explosive  cube  during  cutting.  The  cutter  was  a stainless  steel. 
Serrated  edge,  meat-cutting  band  saw  blade  adapted  to  a stainless 
steel  butcher's  saw  frame.  if  future  applications  and  quantities 
warrant  the  expense,  the  explosive  can  be  cast  directly  to 
sheet  form. 

Preliminary  experiments  indicated  that  the  5000  m/sec  IRECO 
explosive  was  too  powerful  to  use  with  aluminum  plates  less  than 
j/H  inch  thick.  Therefore,  IRECO  agreed  to  produce  a lower  density 
version  of  the  5000  m/sec  explosive.  The  modified  gel  explosive 


has  a detonation  velocity  of  4000  m/sec,  at  a densitv  of  1.2  g/c  m*, 
and  a somewhat  lower  impulse  value.  During  acceptance  trials,  the 


modified  explosive  proved  to  be  suitable  for  the  preset  angle  alu- 
minum parametric  study  and  was  purchased  in  sufficient  quantity  to 
complete  that  task.  It  did  not  prove  to  be  suitable,  however,  for 
the  parallel  plate  geometry. 

In  an  effort  to  meet  the  requirement  for  aluminum  parallel  plate 
welding,  1REC0  developed  a third  generation  of  gel  explosive.  This 
explosive  was  being  field  evaluated  at  the  time  but  was  introduced 
too  late  to  be  considered  for  use  in  the  parallel  steel  plate  study. 

Both  the  high  and  low  density  IRECO  explosives  were  found  to 
possess  some  undesirable  features.  The  tough  water  gel  compositions 
lose  their  moisture  during  exposure  to  high  temperatures  (>90°F)  and 
low  humidity  conditions.  Therefore,  for  prolonged  storage  the  IRECO 
explosives  must  be  sealed  in  moisture-tight  envelopes.  In  addition, 
the  explosive  is  corrosive  to  both  steel  and  aluminum  and  cannot 
remain  in  contact  with  either  metal  for  an  extended  period  without  a 
corrosion-resistant  intermediate  barrier. 

At  about  the  time  the  IRECO  explosives  were  proving  to  be  some- 
what unsuitable,  it  was  learned  that  Thiokol  Chemical  Company  had 
developed  an  experimental  medium  detonation  velocity  sheet  explosive 
Discussions  with  Thiokol  and  a visit  to  their  plant  at  Brigham  City, 
Utah,  to  evaluate  some  explosive  samples  clearly  indicated  the 
suitability  of  the  explosive  for  welding.  It  was  demonstrated  that 
Thiokol  can  produce  flexible,  rubber-like  explosive  compositions  in 
sheet  form,  having  detonation  velocities  ranging  from  3800  m/sec  to 
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m/sec.  Thiokol  provided,  at  no  cost,  HO  pounds  of  tin  1800 
m sheet  explosive  for  parallel  steel-steel  parametric  welds. 

1 hi  explosive  sheets  are  easily  cut  to  the  desired  size  with  a knit  e. 
Moreover,  the.  explosive  is  stable,  waterproof,  flexible  ind  not 
corrosive  on  contact  to  steel  or  aluminum. 

1 Eighteen  ;teel  to  steel  welds  were  made  using  1/4  x 5 x 9 inch 
specimens,  the  parallel  welding  configuration  and  the  experimental 
Thiokol  explosive.  The  welding  charae t er  1st ics  of  the  explosive 
were  very  good,  and  because  of  the  ease  of  handling,  welding  set-up 
time  was  considerably  shortened. 

An  additional  170  pounds  of  explosive  was  required  to  complete 
tin  parallel  steel-steel  parametric  investigation.  However,  Thiokol 
priced  the  additional  explosive  at  approximately  $30/lb,  a price  far 
in  excess  of  the  funding  limitations  of  the  program.  The  Thiokol 
explosive,  therefore,  had  to  be  dropped  from  consideration. 

By  replacing  the  binder  agent  in  their  sheet  explosive  witii  a 
less  expensive  material,  Thiokol  eventually  was  able  to  reduce  tile 
price  of  the  explosive  to  approximately  $10/lb  in  small  quantities 
and  to  reduce  significantly  the  delivery  time.  Under  such  circum- 
stances, use  of  the  Thiokol  sheet  explosive  in  parallel  stand-off 
steel-steel  welding  appears  attractive.  The  modified  Thiokol  sheet 
explosive  reportedly  has  a detonation  velocity  of  approximately 
4100  m/sec  and  a density  of  1.5  g/cc. 

The  Thiokol  explosive  is  potentially  the  most  suitable,  if  not 
the  only  low  velocity  sheet  explosive  ever  developed  especially  lor 
explosion  welding.  If  the  explosive  achieves  its  potential  and 
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the  cost  factors  are  favorable,  the  Thlokol  compositions  could 


make  explosive  impulse  welding  more  generally  competitive  for 
field  joining  of  special  configurations.  Thiokol  has  estimated  a 
price  of  from  $2  to  $5  per  pound  when  quantities  in  excess  of 
2000  lbs  are  required.  Continued  development  of  explosives  should 
be  encouraged. 

While  proof  testing  of  the  new  Thiokol  and  IRECO  explosives 
continued,  a granular  explosive  composition  produced  by  the  Trojan 
Powder  Company  was  selected  for  use  in  the  parallel  steel-steel 
welding  study.  The  Trojan  explosive,  designated  SWP-5,  has  a 
detonation  velocity  of  approximately  3600  m/sec  at  a density  of 
1.0  g/cc  and  a relatively  high  impulse  value.  Eighteen  steel  to 
steel  welds  using  1/4  x 5 x 9 incli  specimens  were  made  using  the 
parallel  configuration  to  prove  the  welding  characteristics  of  the 
Trojan  explosive.  SWP-5  is  essentially  80/20  Amatol  (ammonium 
nitrate/TNT)  made  in  such  a way  that  the  TNT  coats  the  ammonium 
nitrate,  thus  giving  some  protection  against  moisture  pick-up  by  the 
ammonium  nitrate.  As  delivered,  the  explosive  is  dry,  free-flowing, 
and  iias  good  packing  characteristics.  At  a cost  of  $l/lb  the  explo- 
sive should  find  continued  use  for  many  shop  welding  applications 
requiring  large  amounts  of  explosive.  However,  for  field  applica- 
tions where  control  over  storage  and  hand! ing  over  extended  periods 
of  time  can  be  expected  to  increase  the  moisture  content  and  related 
problems,  the  explosive  is  less  desirable  than  the  sheet  being 
developed  by  Thiokol. 
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APPENDIX  E 

WELD2  COMPUTER  PROGRAM 


O 


PRV,  -A*  r L 0 ? (TNPUT,TUrPJT,TaP^^=INPUT,r(\PE^  = r'UTPUT,Tfi^El,TAPt.2, 
1 TAP II, TAPI  A, TAPE  11, T APE 1J, ? AP£ 1 3, T AP 12 b, TAPE 2 7 , T APf  28,  T APE2'  > 

C 

C FINITE  ELEMENT  DYNAMIC,  INELASTIC  AXI3YMMRMIC  CODE 

C FOR  EXPLOSIVE  IMPULSE  WELDING  analysis 
C 

COMMON  H TO  < 1 2 > , 1FL  AG, NUMNP, NUMEL  ,NUMMAT , NEND , IP  A NO , NNN  ,NT,TT,C:iT 

1 UPPI  ,mT  ,f.P  , NUMPC,NST  ,RA  , Al.F  A , JET  A,  TNTSEC  , ‘1  AX  SEC  , MS  TOP  , MTSTOF  ,1  I K 

2 Ml  , «?,  M.-»  , *N,ME  ,l  a , ri,  n,  L UNI,  LUN2,  L'JN  < , LUM  4 , L U F S , L 1 JNG  , LUN7  , N T AF 

3 MAG  * N AGOL  0 »LUN8,<CM 

0 T Mr  N S I ON  A J A ( 0)  , A ( ACT  C C ) 

) AT  A < A JA  ( I ) , I = 1 , 9 ) /oH  A G 0 ,6HA  9 I ,ShA  N - , GHJ  A c ,FHO  R E , 
1-GH"  r.  A,6M  S S C,fH  c I A,6H  T E S/ 

DATA  TSTART  , MAIL  , 1KlO</6HSTART  ,3,0/ 

CALL  <Ln<  ( IKI  OK , 3 , 2 , I PA IL , SEC  , I SEC > 

--  If-  31  MANSION  OF  A IN  COMMON  IS  CHANGED, 

FOLLOWING  S3TtI.JG  OF  KOMHON  lusr  ALSO  p£  changed  -- 
<ommon- a j d l : 

DO  1 I - 1 , < ° M MO N 
A ( T ) = C . 3 
1 COM  T I N'J  E 
L COM -A  : 

NNN0=  0 

1 N-5 
I 0=  G 

L 

RE  AO  { I N , ICC  3 ) L L'Nl  , LUN2 , LUN  3 , L UN  A , L UN  9 , L ON  6 , LUN  7 »L  UN ** , NT  APt  , 

. MAG,MAG0LD,nWOPD 
RE  AO  (In, 1301)  RUN  , MES  T , I M AG 
IF  (RUN. ED. TST APT)  GO  TO  PC 
IF  < I M A G • N ; . 3 ) MA-J  = MA3OL0 
k E W T M D « A G 0 L 0 
I9AIL- T 

‘■’E  A 0 (MA.OLO)  L COM  , (HE  0 ( I ) , I-  1 , LCOM) 

NNNO  = NNN 
? E R I N D NTAPf 
GO  TO  7 : 

r*  ________ 

C MEAD  an  3 h p I NT  OF  CONTROL  INFORMATION 

C 

5 3 IF  ( ir-.ST  . EO.  : ) GO  To  60 

READ  (IN, ICG?)  INTSEC,MAXSEC,NTSTQF 
IF  (MAG. ME. 3)  MAGOLO=MAG 
G3  READ  (IN,  10  00)  (HE n(I), 1 = 1, 12)  , 

. NUMNF, NUMEL  , NUMM AT, NT , N PRINT , NP  » NUM° C , NST ,ALF  A, 9 ETA ,DFLT ,RA 

REWIND  NT  A r>  E 

WRITE (NT APE ) ( A J A ( I ) ,1=1,9)  , (H£D  ( I)  ,1=1,13)  , 

. M'JMf'R,  NUMEL  , NUMM  AT,  NT  , N PRINT,  NP,NUMPC,NST,  ALFA,  -3tTA,DELT,RA 

IF ( NWORD. EO. 0)  GO  TO  6 ? 

M DAT  A = 3 3*b 7*MUMMAT  ♦G’NUMFL  *7 *NUMPC* 2* NP  (-9*  NIJMNP 
nTE?  T = ( (.WORD*  1 0 00-NDATA  ) / ( 8*NUMEL»R*NUMNP)  - 1 
IF(  F' TS1  OP. L£.N  TEST. AND. NT, LE»N  TEST)  30  TO  G? 

WRITE ( I 0, 1G 04)  NTSTOp, N T EST 
NTST  3P  = mTEST 
G?  CONTINUE 

MTST QP=NTST  OP 


■ 


t AO  AOD  P^IOT  OF"  CATA 


I PL  AG  = 1 

7C  M "’=  1 ♦ / * O’JP‘1  t.  T 
0 * = N?  *n  :♦  NUP^A  T 
4 4 = 0 * ♦ OLIO^C 
05=04  iOIIOPC 
Nt>=05  Z 

9 7 = 05  » OUP°C 

• 5 ^ = tj 7 ♦ . jmpc 

0 )z  f.  * '•  (J *-1 

j o11  r. 

• • 1 1 = n i : ♦ ? * n p 
Oil  Ml  = 011 -1 
<Ov<=(<0WM00-fill  + l)  / 2 

1 F (I0AIL .EG.3)  G 0 TO  7b 

■A  1 = 2*  KOI/  3 3 

11 -MI  N ; (VI,  N'JMNP) 

M f : = : i * M l 

12= (<ni- io> /22 

1 ? = p 1 NQ  ( 12,0  'PfL) 

7 b Ni.’  = Ml*.11 
0 1 1 = 01  ♦ 11 
01^  = P1  •■♦•'1 
015  = 0 1 1 

o;5=oi5*',*oi 

O 1 1 = ’ 1 1 7 ♦ 2 * P 1 
0 1 1 = M 1 H ♦ • , * M ? 

I 23  = Mi  ! ♦ 1 fj  * M ^ 

U'c  1=02  3 *0*0? 
n ? 1 A - 0 J 1 ♦ 7 » P ? 

I r < I b A I L . F . . 3 ) GO  TO  9 0 

t r (ma ,l  r . hu pup .op. p2.lt .numel)  ifl ag=2 

IP  ( Pi  .L  T . no  MOP  . 0°.  M2.L  T .0  JMFl)  W PI  T E 1 1 0 , 1 J 3 5 ) ON,  OF 

CALL  9 AT  A 10 ( A < 1 ) , A ( 02)  , A (93 ) , A (04)  , A (05 ) , A (N6)  , A (07  > , A ( Of ) , A ( 09 ) , 

1 A ( 41  j)  * A ( 0 1 1)  | A ( 0 1 2 ) , A ( 0 1 3 ) , A(  (414)  , A (01  5 ) , A (015  ) , A (01.7  ) , A (til  1 ) , 

2 A (013  ) , A(  (419)  , A (02  1 ) , A (021)  »K0MM0N) 


9 j 0E')=;.  *',UMMP 

IF  (IFL AG. EO. 2)  GO  lU  10  3 
IF  (I0AIL.E0.3)  OEOO=NEJ 

022  = 021  A ♦ OF  0 * '1  9 AO 0 

023  = 022  *OEU 

024  = 0 ? T 

N 25  = ^2  •’ 

025  = o?  » 

(427  = 021  A 
GO  TO  2 33 

1 3 C IF  (19AIL.E0.C)  OfGP=(<OMMON-911+i>  / ( 2*  M 9 A 00  «■•♦  ) 

IF  (PI  .GE.MJPOP  . AO  J.M2.GE.0UMEL  ) WRITE(IO,130b>  N E <2  3 , O9A0C 

0 2 1 A = M 1 1 

02?=  01 1 *<0M 

02  3 = 0 22 10  E 9 

024  = 02  7 ♦0UM0° 

IF  (IPAIL.E).C)  p 3 = (<0M- 3*  OiJMnIP ) / 72 
0 25  = 024  ♦ 5 4 * P 3 

E-2 


N 2b  = N ? 1 A ♦ N £ Cl  U*  < M3  A N 0 ♦ 1 ) 
f J «?  7 = K 27 

NMA  X - 1 1 » i J 3 T ♦ *4UMNm 

IP  <N*  AX  . ,T  . KOM  ) CALL  £ X I T 

IP  ( M£  .1  T . 1 . 0->.  NEQ9  .LT.  1 ) CALL  £ X I T 

Z j ; CALL  SOLVE  ( A ( 1 ) , A ( N 2 ) , A ( N 3 ) ,A(N4),A(N5)  , A ( Nb  ) , A ( K 7 > , A ( M H I , A ( N 1 ) . 

1 4 < N ’ j ) , A ( , J 1 ) , A ( N 1 2 > , A < N 1 j ) , A < N 1 4 > , A ( N 1 5 ) , A(  N 1 b > , A < N 1 7 ) , A < N 1 4 ) , 

2 A ( N 1 " ) , A ( N 1 •)  , A <N2  C ) . A ( N2  1 > , A ( N2  1A  ) , A ( N 2?)  , A ('4?  3)  , A (N2  4 ) , A ( N 2 c,  > , 

3 A ( 4?  ) , A(N?7>  , I lAIl,LCOMtNll*ufI->£ST,I<LOK) 


m ° 1 N T 0 ' T - I-  ) A T A A NO  RESULTS 


C 

MT  A=  NNN  -i  -NKNC 

'.ALL  J.JT  JtjT  ( » ( 1 ) , A ( N2)  ,.'.<*43)  , A ( N4)  , A (NS  ) , A (N6)  , A (N7)  , A ( NP)  , A ( K-  9 ) , 

1 A 04  1 J > , A < N1  II  , A (•:  2.  2)  , A ( Ni  » ) , A ( Ml*)  , A ( N 1 5 ) , A (Nib  ) , A ( N 1 7 ) , A (U  1 4)  , 

2 A ( *1 1 1 ) . A ( M 1 ) , A (N2C  > , A ( M21  ) , M A,  NMNO  , A J A ) 

C 


1 j : . PC'3*'  AT  ( 12Ab/  41  5 , 4F  10.  J ) 

no;  FORMAT  (AC,  L!l) 

J?.  P 0 3 M A T l 5lr>  ) 

no?  r at  ( 1 ">  1 s > 

1 ; 3 FORMAT  ( 0 ^ H N T A 0 C HAS  INSUFFICIENT  LENGTH, 

. 1H  p>  _ Sr  T AT  , IS) 

1j3;  FORMAT  ( <♦  I M N01Al  POINT  IP  ELEMENT  3ATA  OUT 
, ?h-‘  NCJAL  OATA  IIUCK  LENGTH,  I 1 j , 

. /'Cm  -'LEMENT  LATA  3L0CK  LENGTH,  110) 

KCb  P0RMAT(17H  COOL  OUT  OF  CCHE/,Z1H  STIFFNESS 
. II') 


ESTAPT  CYCLE  NUMBER, 15, 
OP  COPE/ 


HLOCK  SIZE,  I 1C  , 7h  r.y  , 


SC  'POUT  I'lE  CCUP'JT  (COEFF  jSTPINC, STRESS, ST  PSTN) 
II'TIM  ON  c 0 Tf  (4,  U)  , STRING  (U)  , STPESS(  4)  ,STPSIN  < 4) 
3 5 0 P v : 1-1,4 

0 u h - : . 

00  4 . J = 1 , 4 

'!|H  = CO.FF  (I  , ) ) *STR  i NC<  J ) OUM 

♦ : contin-j11 

c"’M'-iM(i)-nLM 
5.  C ON  T I MU  r 

oof:  1*1,4 

ST  >SSS ( I)  -ST SS(I ) +STRS  IN  ( I) 

b 1 CONTINU- 

1 P T 1 J R ’ ! 

:n  ) 


E-i 


T 


JIPnuT  IUE  CON^CT  ( S T-'flIN  , STRESS,  PROP,  COEFF  , I X,  M , DEL  T , 

. I’J’-lvi  ;f  y l S3  A T ,N:L  ) 

1 1 ME  NS  I CM  :TFAIM1C),STPESS(5)  , PROP  (7>  ,C0EPF  (^>3)  ,IX  <E»>  , Y L J D A T ( , 1 > 
OOMMUU/L  : ARC./LM  ( H ) , .,S  < 4 , <D  , XC  , Y C , S < S , 8 > , C ( 4 , 4 ) , K'  °S  (4  ) 

3 ;N/'L°!.S/''LKfSHP  ,C1jC2»C3,  0 JM, UMAX, &M AX, YCOF  < 2C ) , IIP. CP  ( ! C ) , 

. OOF  ( 1 ),:LbTyti(4)  ,StRINC(4) 

T0  CALL  E L c L 

do  i . : r = 1 , 2 : 
ycqf ; i > =co  :rr  m 
oo  i.i  7 - 1 , s : 

< 0 1 F ( I ) = 0 0 r F f (!♦?:> 

^ 1 J T : 1 ) i . 

,COF  < I)  -'Jjl  FF(  I *50) 

*L  P ( t ) 

v npr  r-pn  P ( 2) 

~.i  = P ■' J-  M ) 

0?=p  J0  Mo) 

»r  PP  1-  ( 7 ) 

' 11"  1=1,4 
T 7 I NO ( I > = 0 EPS  < I ) 

_ ST  ( I)  =CT  * A i *!  < I +4) 


I ,Ml  - 3. 

"JM  = V’AIMFAINI  ))  ) 

UMAX  ■ A 1 S ( C T 4 A I rj  ( 1 1 ) ) 
fi^X^ST1*  SS  (’3) 

Ai_L  tlPL (STRESS, STPA I N,C, SIGNAL, FACTOR, 0ELT,1EPS) 

I F (SIGNAL  .LE.C.)  GO  TO  123 
U'JMYL  0=  NUMYL  ")♦  1 
YL  10  AT  ( 1,  J 0*1  YL  0 ) = N E L 
y L J3 A T t T, NJMYL D ) =F  ACT  03 
123  3 0*4  T I NU  = 

I F (0UW  . AT  .UMAX  ) ST  £.  AIN  ( 1C  ) = STRAIN!  9) 

STRAIN! 0)  = ( 1.*ElSTPN(1)  +cLSTRN (2) ♦ El STR  N ( 1)  *FL  S T RN  (2) - • 2 E * 
. LST-N(4)*  LSTf N (4 )) * ( 1 . ♦ EL STRN  ( ?) ) - 1 . 

TRESS  ( S)  =°MAX 
1-3  3 CONTI  MJE 

30  ?„  ? 1=1,4 

23  3 CTRATNI  !♦’. ) = -’LSTRfM  I) 

RETURN 
F MO 


1 JO 

i : i 

« ■> 

■i.  J ~ 


COPPQ (» , 7, COTE , XO,  ) , X 1 , X2  , I » J , < , L , NO p M , NL A S T , M 1 , NNPT , 


IiJ^ROUT  IN  E 
1 1 '•! ) 

■111  MSI  CM  ^ ( 1)  ,Z(l)»CODE(l),M(l),XOa),Xl(l)tX2(l> 

J ■IMOrvcnRD/  >1  , K J , RK  ,RL  , ZI  , Z J,  ZK  , Zl  ,00  J ( 4)  , 0TS°  ( ) ) , V F L ( 1 1 , A I'iLlM 
T I = I 
J J-  J 
•v  / r K 
lL-L 

1.  |MI  ►•=  NORMT  * M i- M !♦  1 
‘ . M A X - ^1  *'4CPNT 

T F ( ; I . 3 . J ) ,0  TO  ?C 

IF  ( ( r T.  Nxi.j)  .OR.  ( I .ST  .N  TAX)  ) GO  TO  20 
I-P(  I - N M I N ♦ 1 ) 

' 1 - Z < I - NM  I I + 1 ) 

ookd  • . : o o < i - n m i -j  + 1 ) 

.N-r»  ( i -'i'll n *i ) 

1 1 G c ( 1 » = XO ( NM-  ; ) 

JlGD  (?)  = X 0 ( NN) 

/ - L ( 1 ) - X 1 ( M N - 1 ) 
y/EL  ( ?)  - <1  (ON) 

• < COL  (1)  -X  ? ( NM-  1 ) 

AUCL ( 2 ) = X ? ( NM) 

1 1 = C 

2,  IF  (JO.  0 . 0 > V TP  7 0 

I F ( ( J .IT.  n 'IIM)  . QP.  ( J.GT.NMAX)  > GO  TO  n 
J- P < J - N1IM  + 1 ) 

L j=z  ( j-N^iri  + i) 

)O0  ( 2 ) - on  0r  ( J-NMIN  ♦ 1 ) 

IN-  2*  ( J-VTIN*  1 ) 

0 I S P ( 3 ) = XO ( NM-  1 ) 

I I '>r  ( 4 ) =<0(  NN) 
i/rL  ( 3)  - <1  ( IN- 1 ) 

J l L ( •.  ) - < 1 ( N N ) 

301  (3)  - < 2 ( NM  - 1 ) 

ACCL  (4)  =x?  ( N7< ) 

1 J-  0 

3:  IF  ( K < . ).  •.  ) ,u  TO  4 3 

I r < k.lt.ttin)  .or.  (k.gt.nmaxd  go  to  40 

, < = M <-  N'l  IM  ♦ 1) 

Z<=Z ( <-N  1 1 M ♦ 1 ) 

30(3)-  Of)  0:  ( < - MM T N + 1) 

IN  = ? * {<  -MM  I r.  ♦ 1 ) 

I 5P ( 5)  -X  O ( NM-  1 ) 

1 1 j P ( F ) =XO ( MM) 

V-.L  (O)  -XI  (NM-l  ) 

>/EL(G)  = xi(nn) 

A COL ( O!  -X  2 ( NM-  1 ) 

AOCl ( 5 ) -t 2 < NM) 

<<-  0 

4 j IF(LL.EO.O)  GO  TO  GO 

TF(  (L  .’-T.NMItl)  .OR.  (L.GT.NMAX))  GO  TO  53 
•'L  = P (l  -NM  I M ♦ 1) 

ZL  = 7 ( L - NM  IM  ♦ 1 ) 

COO ( 4 ) - CO 0 1 ( L - NMIN ♦ 1> 

N N = 2 * (L-HIINtl ) 

0ISP ( 7) - > 0(NM- 1) 

OISP ( *1 -<0 ( MM) 
j :i I?  >-xi inn-: ) 

E-5 


7CL  ( H * - <1  ( I'  ) 

'\CCL  ( 7 ) = * l < f'*»-  1 ) 
iOOL ( H)  -* 2 CNN) 

L L - C 

5C  'iCHEOK^II  ♦ JJ*KK*LL 

TF(KCH-:  0.  j ) PET  UPN 

If  (U"P‘;T.LT.NL4ST)  GO  TO  GC 
' i * J N I \'N  ° T 

iOPf  T = 

b:  'iC^t'Ts'iCR'J'1’ ♦! 

•£  a r;  ( o \ JT ) c-’iiii)  , 1 1 1 = l , MM » 

',o  m i 3 

'.N3 


E-6 


U 'H  ' !J  T TNE  DAT  A If.  ( P-  >P,  CO  OF  ,H  I , HJ  , V I , J ) , T , 1 hi  , Jf.  J ,P  ,o  , Z , CO  - , : , 

. *0,  X 1 , < 2,  I X , -.1  , EOC  , SIG  , YL  DO  AT,  KO'1vIOf.'> 

•,0‘1F  UN  Hf  D ( I , ir  L AG  ,NUMNc  , fJUMEL  ,NlJMMAT  , ME  DO  , T9  A t 0,  rtfiN , N T , T T , OF  l T 
1 V’P  I NT  ,NP  , NO* PC,  1ST  ,PA,  ALF  A,  Tt  TA,  I NT  SIC  , T A /SEC  , t P',  T OP  , w 7 7 T - , I I k 
? "n  , "2  , * < , Mf. , M£  ,1  A , IN,  10  , L'JNl,  L UN2  , IUN7  , L U fl  U ,LUNb,LUf.'b  , L0N7  , NT  AFC 

j * A 1, , M A f,  ! L 1 , L U N fi  , K 0 M 

I IF-'S  In‘l  PROP  (7,1)  j CO  £ F F ( bL  , 1 ) , Hi  ( 1 ) , H J ( 1 ) , VI  ( 1 ) , 1 J ( 1 ) , T ( < ) , 

. ; MI  < 1 > , J-:  .i  ( 1 > , P ( 2 , 1 > , m > , z < 1 > , C 03 £ ( 1 > , F ( 1 } , y 'j  < l ) , * 1 ( i > , / > < i > , 

. I x ( b , 1 ) , PC  ( 1 . , 1 ) , s I .(0,1)  , y L J 0 A T ( z , 1 ) , P D ( 1 ) , I X F ( 1 C ) , I x N ( > ) 

! F ( ITL AO. E j , 1)  GO  TO  10 
£ W 1 N 3 L'JN? 

£n  I M j L'J  N T 
i:  1 ON  T ! NU r 


AD  A ' . 0 Jj1NT  (;  F OAT:  ' I A L PPQP-PTIES 


0 b t •*  - 1 , IJU1M  AT 

£ AO  (b,10:i)  MTYPE , (PROP ( I ,MT YPE)  , 1 = 1 , 1 ) 

IMP-  OP  <■'■  , NT  YP:  ) .GT  . . . ) R A 0(5  , 10  Jb)  ( C 0 £ F F ( I , M T Y& - ) , I - 1 , c 0 ) 

TF(P -OP  (G,mTyPE) . ,T  .u.)  F F A D ( 5 , 1005)  (CO£FF(T,PTYPF)  ,1  = 21,  ') 

IF  ( F '0°  ( ! ,HTYPP)  .GT  . £.  ) Pc  AO  (b  , 15QG)  ( CO  = F F ( I , h T YP  ' ) , I - ' 1 , ‘,C  ) 

0 i i ON T I MU r 

IT . ( NT  APF ) (( PR C P ( I , M ) , 1-1,7)  ,M=1,NUMMAT)  ,((CO*FF(I,“)  , T - 1 , b 0 ) 

: , m= i , nummat ) 


1 c A 0 A*  J ° IMT  UF  FOCAL  POINT  DATA 


J°nw  - Z 

6 0 P r AD  ( 5 , 1 j j ? ) N ,CODFN,PN,ZN,  XORN , XOZN , X 1 PN  , X17N 
* :jr  F'N  ♦ -’A 

IF (L  . r 3 . Z ) GO  TO  85 
7 X = N - L 

D D - ( PN-PL > / ZX 
07= ( ZN-/L) XZX 
8b  M L = L ♦ 1 
7Z  l=L+1 

I F C N - u ) 1 J 0 , 9 0 , 8 0 

8 j ; o o ► l - : . 

L = P L * )P 
ZL  = ZL  ♦ 1/ 
x ORl  = 0 . 

XO?l  = j . 

X 1»L= C . 

Xl7L-0. 

GO  TO  *1 
9T  PL-FN 
7 L - Z N 

OOOFL  = code.n 

XOPL - XOPN 
X07l - xO/N 
X lRL  r X 1 PN 

x 1ZL  = X l zu 

91  flPOws  MROW+ 1 
I N- N N ♦ 1 
( NF OW) =PL 


E-7 


Z <)F  'H)  -ZL 

GO  ( n • ow  ) ^ r ODE  L 

<n  ! ?*nvow-i  ) = xofl 

<0(  2*N-‘GW)  = V07L 
X 1 ( 2*  N •’OW-l  > - v 1PL 
X1(2*N*QW)=X1ZL 

j.vi)  go  to  tz 
I r ( l . l t . N > go  to  r c 
I r < , -jr>_N»  1 o'„,o?»->; 

13  3 WRITE  (3,20  09)  N 
;sli.  r i I t 

i r F 1 1 f L A G.  iQ  . 1 > GO  T 0 1 10 

- : t ; j i ) { ■ ( i ) , r- i,mn) 

N»0>=  1 

I r ( L . L T . I > TO  7 C 
I » C.UM  JP-N)  i -.3  , 121  , (>: 

110  ;noTit.J-: 

12:  FONT  I N JE 


FAD  OO  : »=>INT  OF  EL  c NENT  NODES 


N = 3 

M3AND= 

' J ® ' 1 In  - C 

1 * ; jE.vi(b»  1 : i 3 > u , ( : x m < 1 > , : - 1 , 1 0 1 


TET 

E 3 MI 

¥ ■ r 

HAN  ) W ID  T m 

I F ( 

ix  ri( 

•>) 

.OF  . Z 

) GO  TO  1 

30 

00 

1,0 

1 = 

1 , 4 

DO 

IGU 

j= 

1,4 

I ABC 

(T 

X M ( I ) 

- I X H < J ) ) 

IF  ( 

'1M.  > 

T 

• 

•V  ) 3H 

- M'1 

COOT  I NO 

E 

L.0 

TO  1 

'7 

0 0 

1 7 f 

I = 

1»  1C 

1 F ( 

( I . 

• • 

3)  .00 

. < I . t U.fa  ) 

) GO 

TO 

133 

00 

1 ir 

J = 

I,lr. 

IF  ( 

( j . 

i • 

3)  .OK 

. (J.ri.n) 

) GO 

TO 

1 32 

tabs 

(I 

X M < I ) 

- I y H ( J ) ) 

it  ( 

T. 

3 1)83 

; M M 

OP) 

TI  NO 

r 

13  3 CONTINUE 

137  MO=?*M‘,+2 

IF(MO.GT.MRflNO)  H 0 A N 0 = M T 
14"  N = 'l  + 1 

TF(M.E  J.N)  GO  TO  14  3 
00  142  1=1,4 
142  I XN (I  ) = IXN(  I)  ♦ 1 
GO  TO  1 40 
14-3  )0  14  3 T = 1 , 3 

148  IXN ( I ) = IXM( I ) 

140  OCNTTN JE 

MKOwrN^Ow*  1 
DO  132  1 = 1,  r 
132  I X < I , NROW  > = 1 XN ( I ) 

IF (GROW ,£Q.M2)  GO  TO  133 

K-8 


o o o 


IF  ( n.  r ; .n  J IL'L)  00  T u 7 0 J 
1 F ( >. . j .y ) r,  ) TO  i .• : 

. ' TO  i >•. ; 

! • • • : F ( I f t : , . • , . 1 ) .0  TO  ' 0 1 
•J  T j:>  ' ) ( f 1 > , 1 - I , 'IE  ) 

4 'IT-  ' ■ : T A F>  ) ( ( ’ ' ( J , <)  , J = 1 , r-  > , < = 1 , 1?  ) 

,>'h:  : 

I - ( f . . . N - : L ) V)  T " r ] 1 

i r t r. . . . m O ’ t : 7 , 

t j i : 

r . : . nT  I \ i E 

.i)  rp  i t - (Luoi)  (E n ( i ) , i - 1 , ,<r i 

: * ( • r ■.  -»  > < ( i * ( j , o , j - i , o > ,<  = i , > 

L . J-  > . .,V.  A o y OQ.'I  ) I T IONS 

7 0 E . r O.  lie  . ..  j)  J TO  4 J 0 

ETC  < - 1 , t !JMO' 

* •:  A (‘.,13  7 ) T N I ( K ) , J N H < ) , A , H , T ( < ) , ->  1 , P J , Z I , 7 J 
I = INI  ( ->  ) 

J - J N ) ( < ) 

0 z - ( z i - / j > / 1 . : 

0’’  = I*  j-p  i > / u.  c 
■ x - a * ( 7 . : * 7 : ♦ ; ■ j > ♦ ) * ( - 1 * » j > 
z * 

Mi (k>  = = x*  DZ 
H J { <)  - 7 < * 0 7 
y/I  (M  D ~ 

1 J ( K ) :/X*  jP 
7 ' . ONTINU' 

^ - 1 T : ( NT  a PE  > ( HI  ( < ) , V I ( K > , H J (<)  , 7 J ( K ) , T ( « ) , IN  [ ( K ) , J")  t ( k ) , 

1<= 1 , NUMCD) 

P'A'j  Ar.D  p P If  T OF  LOAD  DATA 
) 0 / w G w = 1 1 N P 

7 . PEAOC5,in«.j  ( P { < , M ) (<  = 1 1 ?) 

W‘-IT  ( (TAP  ■'  ) ( (P  (K  , M>  , <-  1 , 2)  , M = l , N») 

<♦0  0 CONTINUE 

IF  ( IFlAG.  ► J .2)  PETUPtJ 

NO<-  G Z*NUMM  ATf  1 4#NL'NNPf  2 2*rJUMEL  ♦ Z*NP»7*N'JHPCf?*  NOMNP*Mn  AND 
I F'  ( N 0 . L .ROMMON)  F E TU°N 

1 F L A r>  = ? 

r ^ I ' 1 3 INN? 

• £ W I m 3 LUN? 

^P I T - ( L UU  D ) ( P ( I ) , I -1 ,Mf!) 

W-  I T r (l  J N 3 ) ( . J ( I ) , I =1,M£) 


RETURN 

1331  FORMAT ( 13 ,4F15  . u,  3 F *5 . 0 ) 

10CZ  FORMAT  (15, F5.  1,6F1C.0) 

10  2 3 FORMAT  ( 1 1 I E3  > 

IOC*  FORMAT  ( .?  F 1 C • 0 ) 

1C  05  FORMAT ( 1 ] F S . 0 ) 

100  7 FORMAT  (?  15 , 7F  u.  3 > 

?0CD  FQsmjt  (DOHONOAL  POINT  CARO  F 7KOtf  N=  I r> ) 
END 


!»•--  t L-  CT  ( '1  , r, , c > 

jime’JSIO'i  r,  (l,  ui 
C ( 1 , i ) - T ♦ 1 . ? 3 3 * r. 
CIl,’)  = \-.hbb7*& 

C ( 4 , 4 ) = n 

r,  ( i , u - : < i ,2) 

C ( 1 , 4 ) - o . 

C(2.1)=C<1»2> 

C<?,?>  =C  < 1 , 1 > 

C(2, i)=C(1.2) 

C ( > , 4 ) = 0 . 

CM,1)=C(1,?> 

C(T»2)=C(1|2) 

cn,5)=c(i»i) 

C < 3 » 4 ) = 3 . 

C (4 , 1 > = 0. 

C < 4 , 2 ) = 0 . 

0 ( 4, 3 ) = j. 

UPTURN 

ENT 


K- 10 


SUBROUTINE  ELF  UN (STRESS, STRAIN) 

JlMENSICN  STRESS(l)  , STRAIN! 1) 

• UMMNN/ El PLS/S , G,C  1 , 02 , 03 , JUM, DM  AX , PMAX , rCOF  <2  J > , 9COF ( T 0 > , 

. . 0~r  < 1 ) , EL  r • \ (4  , , s TRINE  ( /,  ) 

I F (~  : .*_  . „ . ) GO  TO  1 J ] 

( VAR!  V‘Lr  A < Mi  Gil  IS 

lsnjnci 

_ 2 = • _ ' T R'  I { ) 

< - r . c ; . j , j j 

* 4 ■ E - ' T R * , ( 4 ) 

, : - r.  r - \ ■;((:.*  : + 2 + E l » - z - , 2G  * L 4 * c 4 > * ( i . » e ; ) - i . > 

: - 3 0''F  ( 1 ) ♦«.  . » »(C°F  ( i.  ) * JCRNT  4-  3.*  3C0F  ( 3 > * UCRN  T*  UO  P N T 
1 Ju  IF  (CT  ,L  - . 3.  ) RETURN 
C VARIABLE  SMEAR  MODULUS  • 

BUM? = STRESS  ( 1 > +STRFSS(  2)  ♦STRESS  ( 3) 

IF(DJM?.LE.PMAX)  GO  TO  2C0  ; 

C UnLOADING/RELOADING  j 

G = 3 C 0 F ( S ) ♦0UM2*SC0F  (6) 

IF (G.GT  .SCOF  ( 7 ) ) G = SO OF  ( 7 ) 
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. 1 x c , i ) , :p'.,  i i ; , i) , sio< g,  1 , yl  3 dak  2 ,1  > , ed<  i>  , a ja  (O) 

■’Euro  K!  T A •’  F 
■|!p  L K-  (UUUMr  *■  Ml  - i ) / *u 
M El  L u'.  = r ! ijm  Et.  4 M,E-  1 ) / M2 
r e ( . . 0 . ; ) go  n 45 

'■>  C A ...  ( ’ j ' A ) (AJA(I), 1=1,0)  , ( HE  D ( I ) ,1  = 1,12)  , 

. UUM;  5 , VJ'-’CL  OSU-UMAT  ,NT,NPMMT,NO,NUMPC,NST,  ALFA,  BETA,  OEl  T,RA 
4 ct  t ( 1 u , 1 0 ; :■  > < a j a 1 1 > , 1 = 1, 9 > 

u=>ite  ;n,co::>  ihekd,  i =1 , 12) , nuhnp,  numel,  uummat  ,nt  ,e  o^int  , 

. ?40,  NU'^’C,  ALF  A , V T A, DEI T , UST 

7 = r.‘  r a 

r F ( UN. JO.  N-  . 0 ) on  T(j  19  ] 


-EAU  AND  P.  I N 7 OE  WA  7 cHI  AL  opqPEPTIES 


1 


it  AD  ( rO  A • : ) ( <P >C°(I,M),1=1,7)  ,M=1,NJMHAT)  , ( (CCErE  ( I , M > , I = 1 , 6 G ) 

1 , M= 1 , AT  ) 

-(PIT  E (IO,75u:>  (ADA  (I),  1 = 1,7),  (HED(J),J=1,12> 

t - < • , , i;/i 
D ■ I M = 1 , ULMU  A T 
1 T Y P " r 4 


UP II  - ( > , ■>  t . c ) MT  VPE  , (R0D  ( I , "I TYPE)  ,T  = 1 , 4 ) 

I F (POO-  r , M T T p E ).>!.:•  ) WRITE  ( G , 213  0 (C0EFF(I,MTYPE),l  = 2i,r'C> 

1 f(p  UP <7 , -TYPE ) .GT.  :. ) WP I T E ( 6 , 22  E C ) ( COE  EE ( I , MT Y P" > , I =r  1 , 6 7 > 

L N-:l  a 3=  Pu0:1  (E> , MTYPE  1 

IE  ( IUEL  A-s.LE.:  ) .0  TO  GO 

;F(Ii:El-VO  (j.l)  WRIT-  (6,2300 

IFdMLAS.  C .?.  ) WK  IT-'  (F  , 24"  Z) 

iE  ( I Mill  A 7 . E 0. 7 ) '-O  IT  i ,25 *3) 

W P I T { b , ' 6 j J ) (Or  EE  (I, MTYPE),  T = l, 20) 

;lut  i nj-je 

4PIT'(',7Gvr)  A J A , h ; j 
wPI  TO  < , , - .11 

)C  i . l~  1 , r:-L  ol< 

•’ETC'  V U T A p • ) ((I  < ( J , < > , J = 1 , 5 ) » < = 1 , M 2 ) 

J S ’ A > T = v 2 * ( I - 1) 

USTOr  = '■  'JM  L - IT-  T .IP  T 

jsr l p m r o ( m ?,  uo.f  e ) 
ir  i,r  K-i,MOTrp 
M E L = < ♦ I . i T A T 

wRT  T • (S,2.:3)  NEL , (IX ( J,<) , J=1 ,5) 

TOUT  t mu 

IFIMU^PO. lU, , ) GO  t0  1*3 

•’■M.  (UT  A'  E)  ( «1  ( O , 31  ( <)  , H J<<>  , V ) (K>  , T ( K>  , I MI  ( K)  , JN  J(  K)  , 
i<  = ; , mump'.  > 

AO  ('.TAPE)  ( (°  (K,  0)  , K=  1 ,2)  , 0=1  , NIP) 

URI  T ■-'(-,,  ) AJA,mEO 

u u T T - (6,  ’310 
10  ' r C <=1,NUMPC 


I- -21 


i-im  ; «. ) 

J - J N J ( < ) 

•••.  4®ITr  ( S , J j l 3 ) I,J,  Hi  ( K > ,VT  <<>  ,HJ(<)  ,v  J ( K)  ,T  (K) 

01A1  W ' ^PIilT  Or  LOAD  '3  A T A 

■I  •?  I T ( - , ' s . J ) A J A , h e o 

w w ■ t r ( 6 , j _ r r > 

T 3 I T T (3,?3G5)  ( (?  ( <,H)  , K=  1 , 2)  , H = 1 ,'JD) 

i h : io  i'll  mk- i , rj r l k 

3 r A ( ■ ! T A o _ ) ( ( s 13(1,0)  ,1  = 1,4)  , ( EPS ( I , N ) ,1  = 1,4)  , N= 1 , *2 ) 

1 H 1 1 •‘ITIi'iiJ- 

T ( , . , ) A j.l  , H£-| 

W'"lTr  ( , 1 3 . 4) 

■>  i 130  \ 3 < = i , f 1 1 , • ilk 
3 E A 0 ( M T APE)  ( 3 ( I ) , I = 1 , ON) 

M iTJ^T- •'1*  ( N OK  - ; ) 

jSTv.-  > = (.-/  3 J*3  -ms  t apt 

MC1  C'-'=  3 IM  ; ( Ml,  '.'STOP  ) 

V3  1 1 j K = 1 , IV  S T 1 ? 

IM  3 = v t"5Tft-’T 

• i - I T ( . , J c)  • ■M'J,COr!E  (<)  , 3 (<)  ,Z  (<)  , X0(2*<-1)  , X0(2*0  , XI  (2*«-i)  , 
. x l ( 2 3 < ) , 1 2 ( ? * < - 1 ) , X 2 < 2 * O 
i . : oo*i  t r nu“ 

1 »;  ;o*.t  iM'.'t- 
T T - t-  f 1 , 1 ) 

- L = G 

I r ( IM  !C.  £ . . J ) -0  Tr  100 

yri  jMr  = "r/jo 
1 T - T T , i ( ,/j  ■ ji , * >0  L T 
_L  - m 10  ( i'.MIO  , MPP  xMT  ) 

: v,  jC  o , j m n:=  l ,mt 

•-  L : L L ♦ 1 
M - T ■ O : !.  T 

:r  ( L L . _-..N°PIMT)  r.u  TO  3 C 3 
in  P . C 1=1,  NflP  31  < 

LAO  ( .TAP”)  l w < O , 7 ( < ) , X 0 ( 2 * < - 1 ) , X 0 ( 2 * < ) , X 1 ( 2 * < - 1 ) , X 1 ( 2 * K ) , 

. x 2 < 2 * < - 1 ) , x 2 ( 2 *<),<=!, Ml) 

,’3  1 'A.  IT  IN’ Jl- 

JC  V • : 1=1,)  7 L 'L  < 

-'S..n  ClTAC-  ) ( r I -J<  J,N)  ,J  = 1 ,4)  , ( Ep  S ( J,M)  ,3=1,0  , N=  1 , M 2 ) 

?S’  :0NTIN'J  = 
r>0  T *< 

3,  j L L = C 

^rit- (- , :i  a ja , he o 
4 p i t (-,?::<>)  tt 

?::g  i = :,n‘iP  )l< 

-'SAL  (.T',r>)  (=»<K),7(K),XO<2*K-l>  , XO  (2*  < ) , < 1 <2*<-l)  , X 1 ( 2 * K ) , 

. < 2 ( 7 * < - 1 ) , X z ( 2 * * ) , K = l , M l ) 

MS  TA^’T-f'l*  (I-n 
M S T C c = ' jMflo-MST  AWT 
MS'  ■ > = 'T  IM  ( NST  OP,  ) 

JO  3jG  <- 1 , NST  0° 

)N°  = < +MST  A-'T 
A A=  x 0 (3*-:-l)  *P  ( K ) 

3 3=  X )<:*<>  f ? ( K ) 

ww[  T-  (6 ,2G3fl)  NMP  , AA  ,33  , XO  (2*<-l ) , XO  <2*0  , xi  ( 2 * K-l)  , XI  < 2*0  , 

. x 2 ( 2 * < - 1 ) , X 2 ( 2 • O 


K-  22 


CrNT  I NNr 

M R L T (.'■,?  5 j C ) ' 

W P I T (:>,12j£)TT 

writ: , i : :o) 

: l-  i,nel  j 

‘J  ST 

p = A D ( J T 1 p r ) ( ( 

\J,'  TP  » ; ‘.J-J  . f 

*!C'C  13  = '1  I ) I < M2,  V 
TO  u j j <-  , r.ST 

j r L - + N ■ f ART 

writ  iio:> 
COM  I N j - 
com  • j : 

F TU  VJ 

F ORM  a T ( 2 J X , ic H 
-0'M,U(n«  CL.NO 
. /*,'.  p -R,77,5 
P OP^  : T ( ^ T I 1l 
F C p « A T ( I . , , * ► 1 ? . 
f-'O-'y  at  ( i h i , ?:x  , 

r r . T l!i(.17/lh 


, ?5jC>  '•  JA,hFO 

,12jOTT 
,1120) 

I - 1 , ML  Cl  k 

M2  * ( I - 1 > 

T A ' ) ( (CIjIJjM)  , J = 1 , 4 > , (EPS  < J ,W)  , J-  1,  -»)  , N=  1,  M2  > 

J'/l  E - N'T  1 T 

I J : t M2  , T)F  ) 

< - , f.ST  m< 

T r A R T 

1 i 0 C ) MEL  , ( S If,  ( J , <)  , J=  1 , 4 ) , ( £P5  ( J , < ) , J-  1 , 4 ) 


TARLE  OF  MATERIAL  PROPERTIES  //) 

. 7X  , 5HSIG-  R,  7X,  5HSIG-Z,  7 X , 5HSIG-T  , bX,  1 0 T AC  - ' L , 
• ■ _p 5 - Z , 7 7 , 5 H pp  S - T ,bX,bHApS-RZ) 


H i TIC  T -E  1 0 . 5 ) 

I H 1 , 2 C 7 1 ! A 0 / ) 
(17,127  -./ 


i J *i 

» 4 m "j  r ^ 

OF 

NOCAl  POINTS 

14 

/ 

o 

V" 

• MIR 

OF 

-L  I H-NTS 

I 4 

/ 

■ u nr  - 

OF 

01PF.  MATrRlALS 

I 4 

/ 

"4 

j-'HE  i< 

n r 

TIFF  INC  P EM  ENT  S 

14 

/ 

C 

7 >-  ~ 

P J I NT 

IN  T 

R V A L 

14 

/ 

1 

< > 

f J M .)E  P 

0 F 

COAT,  POINTS 

14 

/ 

7 

7 2 ^ - 

, J 1 1 •.  k 

Of 

- E STUPE  CARDS 

I 4 

/ 

i>; 

AM  PI  (j 

& 00 

EFFICIENT  ALFA 

F 1C 

• 

5 / 

’ ) 

3 ' H 1 

; AH  PI N 

r r'  -\ 

EFFIPIENT  U ETA 

FI] 

• 

5 / 

t 

T 1*1  I INCP- 

ME  NT 

E 1 2 

• 

5/ 

• 

7 ' M ' 

* ; u or  p 

o r 

VANISHING  EL  EH 

I -j) 

r ORM AT 
FOR* A T ( 
F 0 3 P‘  A T 
F 0 7 N A T { 
, L P- 


T UOH  r L f N r N T NO.  I 

T < 17 ,F 1 J . c , hF  1 C .3, 2ElC  . 2) 
T (111.3,410,  11?) 

T ( 7 M H NO.  V.  POINT  T y o £ 

P-  7'.  L /-7-.L  , 2JH 


M-ORD 
p-  ACC 


Z-ORO 

Z-ACC) 


M A T E 0 I AL 


R-  IT  SL 


2 C 1 3 

->  , i * 


r QPh A r 
FORMAT  ( 
. H * , L m Z- 
FOOf  A T 

F C RMAT  ( 
F 0 P M A T 
1 5 X , l h I 
0 2HV  J , 1 
F C JM  A T ( 
P P ^ A T ( 
FORMAT ( 
F 0 R M A T ( 
F QP*-  AT  ( 
FOR*'  A T ( 
FORMAT  ( 
FORMAT  ( 


(?F 10.7  ) 
(H  T I ME 
1 >°  t )*  > 

( VM 

10 , MF14  . 

( O DP  r 

, '*  x , 1U  , 

1 X , 1 W T ) 

_ I t>  , 2 4 X , 

1 ’■  M MUL 
1 3 M 5 HE  A 
1 7 H 01" 
13H  NEW 
1 7 H VISC 
13H  y IEL 
?4H  m AT  p 
nUL  < 
?4H  SHfA 
24  H QENC 
.h  TmIC< 


T-E15.5/5H  NODI , 9X,5HR-0RD , ) X , 5H  Z -ORD  , 8 X , PH-‘  - 0 I S F , 
OHR- VrL  , )7 ,5  HZ- VEL , 9X , 5H 9- ACC ,9X ,5HZ- ACC) 

TIME  PRESSURE  D) 

5) 

35U  PE  HOUNTARY  CONDITIONS/ 

M , h H PI/P, 3 7 , 4H  FJ/P, 8 7 , 2 H H I ,107 , 2 H V J , lCX^HHl,  1CX  , 
:.F  1 E . 4) 

< COF-FP  ( , II  E 1?.  4 , /l  TX  , 1QE  12 . 4/M  x , l :E12 . 4 , 2H  >/) 

" COEFF ( , 1 JE 12. 4, 2H  ) / ) 

PLASTIC  /) 

PLASTIC  /) 

O-PLASTIC  /) 

D CCEFF ( , l IE 12.4/lOt 12.4 , 2H  )/> 

p I A L NUM9FR  15/ 

MODULUS  Fi5.5,oH  (PSI)/ 

p 1C0ULUS F15.5»6H  (PSI)/ 

I T Y p 15. 5/ 

MPS?  (ONE-!) FI  5. 5, 5M  (IN)/) 


>a*://  i:x,i 2kb//) 


0,,>.  4 T ( :m  , 2 jx  , 
= '!  T 


b‘j  V -tj  r [ > I r Pi.f.  :'T  (F  , •),  ,,C  ) 

) I MF  'IS  i -fl  r ( 4 ) |CU*4l 
f i = r ( i » 

= r ( >) 
f 3=f  i i) 

F -«=»  14) 

C lA<:  S 3 A- A-'l  T 
X = 7*  3 

(-P-fl  ♦ r 2 4 F 7 

n=‘ri*>-l  ♦ f 2*f  L;  ♦ F j*n  ♦ 2.*F4*F4 

jrX  *F  r*  Fr  +2  . *F.  * 5 
T F ( > . ■ , . . . ) '1  = 1. 

A1-X*FF ♦?.*G*F  1 
A ? = x « F c ♦ 2.  *0*  F 2 
A ? r rrc*  F J 

r 

1(1,1  </♦?.*•.)  - < A 1 * A 1 / 2 ) 

Kl,  2)  - x-(A**A2/D) 

,(1,3)  ' X- (A1 * A 3 / C ) 

2(1,4)  - -?.*G*F4*  a i/i 

c 

C < 2 , ) = x ♦ 2 . * •-  ( Ac*A2  / 1) 

2(2,3)  - x-  ( a ' * A3/D  ) 

2 (2,-)  = -';.*'."F4*  A2/rJ 

C 

C(4,j)  = X*2.*r'-(A3*A3/‘)) 

•2  ( T , u ) = -2.*(-*F4*A3/D 
r* 

■J 

2 (4  , 4 ) - 2-4.  *2*  (,*F«4*F4/0 

2 ( 2 . 1 ) = C ( 1 , 2 ) 

>2(?,1)=C(1»3) 

2 ( < , 2 ) 2(2,3) 

C(4,l>  =2(1,4) 

(1(4,2)  = 2 (2,4) 

!(,,  < ) = 2 ( 7 , 4 ) 

2,  (ETli-’N 

-:rn 


F.-24 


n o 


3LHR0UT1UE  OUAO  Ul,»2fR:i,94>Zl,Z<>,Z3,Z4,RM,Z'1>VCl.,'),aK,0S> 

c 

FO’fj  STIFFNESS  M A T R I X JK,  CENT  ROI  QAL  STRECS  MATRIX  QS 

F 0 p A F i 7 E POINT  AXISYMMETPIC  IRON  S 0 U AF  =>  I LA  T E P At  USING 
C A "njc  POINT  INTEGRATION  FORMULA. 

C CONS'  A:  iT  '^rv  STRAIN  INTRO  DUG  F S INCOMPATIBILITY 

} I Mr  \I01  O’!  K (rt  , i)  , US  U,  S ) , 2 (4,  rt)  ,TT  (4)  , CJC  ( 4, 1 J ) , SS (4 > , QQ  < 1 C , 1 C > 
jATA  S'>/  -I.»l.jl»»~l./»TT/-1#>-1, 

C 

do  *•  i- 1,  i: 

)0  b ) - 1,  l'J 

- j j ( i , j)  = : . 

..>1  1 r > « _ . o 

R1  ? = 

, ' j - n _ >■» 

•V  4 = ■ 2- J4 

/l  1 - / 1 - 7 2 

l i — / _ / . 

/ 1 - * / 1 - 7 4 

/ 2 < = / 2 - Z < 

/ ’ -4  - 7 P ~ 7 -4 

» 7 = / ■»  _ / 

7 CL  - 1 ’-*7.2*-*?  4*213 

pm=(P1«p->*'  x »■  0 4 ) / 4 , n 

7v,-(Zl»7',*-2  3*/4)/4.Q  | 

if  r>  ( : , 1 1 . l j.  ; > go  to 

YC=Z?4/  y;,L 

x G = - 13/7  ) L 
x / /VOL 
Y-l-Zl  V 70 L 
X 0 - - « 7 

Y G = - Y H 

Y 7 - “ Y 1 
/ H - - x f 

3 0 . 11=1,4 

S=  »S ( J : ) * J. 5 77  'G02fc91fl962G 
T = tt  ( j j ) * ; . 07/  30rj2b91H9b2G 

/J  =7,L»Gt  (P’A*  712-R12*  7 34>  +T*  (R23*714-Rl4*Z2i) 

XJAC-X  J/H,  . 
jM= 1 , C - S 
c=  1. C *S 
T«= 1 . D-T 
T p= 1 • C ♦ T 
h l- c . : 0*3 M* TM 
H2=  . :g  » jp*  T*« 
h 3- C . 2j*  jp*  f° 

P4=C  . 2C'*SM*  TP 

=>  = Hl*Fl*M2*f<2*H3*R3*H4*R4 

0 l=Ml /> 

G t=h  V 
,4=M4/R 

GC=SM*SP*TM*TT-  /p 
xi=(-P2'oRT4*S*P23*T)/xj 
x ?= ( p l 3-R34*S“-14*T)/xj 
X3=(  R24-pi?*SyR14*T)/xJ 
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▼ 


O U O 


2 . * ( - 

20  < - 
2 •„•*(- 
:-,*(- 
' , ' » ( T 

• IT 


1 2*S-p?3*  T)  /X  J 
34*C-Z?3*T)/XJ 
' 4 * S * Z14*T)/XJ 
1 7 * 0 - / 1 4 * T ) / X J 
1 2 * S ♦ /2  3* T) /X  J 
"N*r^l+TM*02*TP*,:,.3-Tt3*R4) 
T’''*/ifTM*Z2+TP»77-TP*Z4) 
.t;p»p2  + r;D,'<3+SM*P4) 

' ” * 7 1 - SP  * Z2 ♦ SP*  Z 3 f£M  * Z 4 ) 

» *RT  ) / < J A 

: S M * .5  P * Z S - 5 * T M * T P * 7 T ) / x J A f. 


r AO-  x j .1  ; * R 
F ST  IF  FM  SS 


n = r <i,i 
) 2 = 0 t J , ? 
1 Z - r < I , 7 


. * •* 
• \ , t 

» L > 

. 7>  *F 
, •’ ) * F 
,-*>*•' 


DC  ( 1 

, 1 ) 

U*  Y 1 + 14’  tb  + 03*0  1 

j C ( i 

, ■'  ) - 

0 1 * Y 0 ♦ J4  * XG  0 7*0? 

i ; i 1 

, r 1 : 

1*Y3f04*<7+05*G3 

j'.  1 1 

, 7 ) * 

1 * Y 4 f J4  * Xft  *030  4 

, ",  1 1 

, i ) - 

' 1 * Y C * 0 3 * G C 

j C ( ! 

, ) = 

•_*X  1 *04*  Y5 

r.  ( i 

,41- 

...  * X 2 * 04  * Y h 

OC  ( T 

. ) = 

jj*X 7 * 0 4 * Y 7 

) : < j 

, J ) = 

.2*  X4  *04*  Yft 

jr,  ( i 

, i : ) - 

n : * x c 

CO  JT 

1 0 J" 

>o  ; 

3 f = 1 

» * '■> 

")  L = 0 

o(i,t 

) 

^ 3 - 

n f ■>  T 

1 

- — 

» f 1 

J 

, > - , 

: ( x , : 

) 

4=  D 

. ( * ,1 

) 

.4(1 

, I ) - 

11(1,1)  * D 1 » Y 1 + 0 4 * X S F 03* 

n < 3 

, ! ) = ) 

j(7,I)  f01*Y2f04*XF>fQ3*G 2 

J ,(  S 

, I ) = 1 

1(0,1  ) * J 1 * Y 7 f D 4 * x 7 f 0 3 * 0 3 

10(7 

,11-1 

1 < 7,1  ) *01*  Y4*04*XftF03*G4 

10(1 

, I ) = ) 

1(0,1  ) *01*  YC  * 0 3 * C»C 

1 J 5 ? , 1 ) - J J ( 2 , I ) ♦ 02  * X 1 f 0 4 * Y 0 
1 1 ( 4 , 1 ) - J 1 < 4 , I > f02*X2fD4*Y£> 
J'i(F  , I)  - i J(F,I  ) *D2*X?f04*Y7 
J J( ft , I ) - J J ( ft  > I ) *02*  X4f04  » Yft 

i j(  i . , i ) = ')j  a: , i)  ♦ o?*xc 
COOT  ! N'j  E 
CONT  I NOE 

FQVf  STOSS  MAT -lx  QC  AT  CL  I 


MAT-IX  QC  AT  CLMTRIOD  (PM,ZM>  0 F I L E Mi MT 


DO  4 ] I -1  , a 
01  = r ( I , :> 

02=0 ( I , 2) 

13=0  ( I » xj  X (4.Q*pm) 

j4s?  < I , 4) 

'!=(  01*/2h-04*p?4) /VOL 
T2=  < -?1*  T 1304*  Xl  3 > / VOL 
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r 

p 

0 


0 

r 


T?z(_-ij*o2_(+i4*Z24)/VnL 

;<♦=<  q?*  u i-rs<.*  zi3> /you 

ao  ( i , 1 » = a i*  t 1 

ICU  « *>  0 3*T2 
)cu,ca=!^-Ti 
3C  ( I » 7)  - ' ?-T? 

) C ( 1 , 9 > = ** . 1 * 0 J 
)0  ( 1 , <r>  = T i 

acn  . •♦>  = T- 

1CU  , C>  =-T  ' 
iC(l  o>  --T'« 

1C  ( I > 1 , ) - C . 0 

■* , coaTivjt: 

cr*iirr.  modi 

DO  C- . 0 = 1,2 
L ~ 1 3 - *• 
m = l ♦ 1 

JO  r : r - l , L 

C=  )0  ( 1 , ,J!)  t 10  <M  , M) 

00  5 - J-1  ,L 

v . 30  C I , J)  =0  0 11  , J ) -C*!jQ  (M,  J) 

=>FLGO  AT  bTPFOS,  STIrFUFSS  AN)  LOAD 


HOH  COOTTN JF 

) 0 f J J - 1 , 0 

jo  7 : 

JK  < I , J>  = 1Q(  1 , J > 

7':  ) K ( I * 4 , J)  = 3 0 < T *h,J> 

301  = VOL  *••*'1/  2 . 

OFT  i;  r ■ 

~ NO 


I 


1ATRICFS 
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. ih'  ' i’r  - :s T a~  ( hail  ,1. com,  prop  >f  n u , h as?, - , p,nnprl 
1 rj"t  pf  , ;,nstpl-  > 

p f a 1 m;,-; 

COM)  *1  i-  ) ( 1 ) , If  L AG  , IIUMNP  , NU’-IE  L ,N'JPMAT  , NE  JO  , MR  AND  , Nf.’N  , NT  , T T , CEL  T , 

1 NPP  1 NT  , \uM  - .'•,MST,RA,flLrr,'JETA,rNTSEC,rtAXSEC,KTSTCP,MTSTOP,IIK, 

2 Ml,  , •"  , " ,L  A , IN,I0,IUN1,LUN2,IUN3,L  JM4 ,L UN5, L'JNG , LUN7  , NT Af  , 

••  Mil  , * " GDL  , L U r > - , K C M 

: IMMJS  ! . • pp  ip  ( 1)  , MA^S  ( 1 ) , E ( 1)  , p (1)  , rD<  1>  , A ( l) 


c 

c 


U'l-  ’*  .U  <Np 

i k < i«a  !l  . . . n 0 j tp  r- : c 
T Z •*)  that:  On  (.•••:  Tn«.'.,.r  :y:ll 


pE'ir  t ma>, 

I r o l ■ 1 1 1 1 ) .0  to  1 ] c 

PE-JT'iP  L'JtJ ; 

- E ■<:  ' ' o i j ■ r: 

I J L'JM  < 

*irid  c':.  ,)  _co:s  c h c j < x ) ,i  = i,lcom> 

WPI  T (-  • ,)  (••>-(  ■ ( I ) , 1=  1 , rmMi  ) , ("(I)  ,1=  l.NEC)  , < MASS  ( I > , I - 1 , NUMNf  ) 


00 

1 

--  oi/iN 

L < 

i F 

( 

' F L A ; . N ■ . 1 ) 

p-;  a j 

(LUNl) 

( A ( I ) , 

I - 1 * LA) 

wt  [ 

7 

('■'.',)  ( A ( 

1 > ,1  = 1 

» L A ) 

• > n m 

r 

I N J E 

on 

1 

♦ V.  “ - 1 , ( NF  R 

L K 

I F 

( 

F 1.  . 1 ) 

pp  AD 

(LUND) 

(R  ( I ) , 

I = 1 , M N ) 

g p I 

** 

(PAG)  ('  ( 

: ) , 1=1 

, MN) 

1 u : eon 

T 

INF 

00 

1 

• . *:  1 ,0"L  :i 

i f 

( 

• 1 ) 

u:  A 0 

( L U N 3 ) 

< L 0 ( 1 > 

, I = 1 , M E 

R P T 

; p a r, ) ( £ 0 

! I > , I = 

1 , HE) 

10:  con 

▼ 

I N “ 

! f ( : o a 1 l . _ ; . 1 ) 

IF  ( I 3 A I L • • . . n 
if  (i p ■ 1 1. . . . r > 

-r  ^1  no  v « . 

U = MAG'H  I 
‘ ! AGO L 0--  ’* A G 
M AG-  N 
RE  TO0 N 

> - S T A R r 0 Y C L E 

vJ  J f ( IFL  A.» . LC.l  ) ■')  TP  02C 

-ERlNI  l'JNI 
•'ERlNJ  LUU  2 
■F/llNO  L’JflT 

‘,2  REA  . (M1VJL  ))  ( pppp  ( I ) , I = 1 , igllMl ) , (E  ( I ) t I = 1 »N:1>  » 
1 ( M A ^ S ( I ) ,1  = 1,  N'JMN  F- ) 

JO  0 * 3 N=  1 * NOT  ">L  ^ 

••  E A j ( 'A  .OL  0 ) < A (I),  1 = 1,  LA) 

IF  (IFlA,.Nc.1)  WkITE  (LUN1)  (A(I),I=1,LA) 

>4;  CONTINUE 

JO  ; *4=1,  NNP  PI  «. 

PEA  ( " A jOL  D ( P ( I > ,I=1,MN) 


■HP  IT:  ( I 0 , 2 o u 0 > NNN  , MAO 
NT  = NNN 

W a 1 T J (10,2301  ) NNN  , -1  AO 


E-28 


IF  (IFuAG  ,NE.1>  WRITE  ( L UN  2 ) <MI>,I=1,MN> 

So . OONTTNU- 

20  E J.  C U= 1 , NEl  Di_< 

REAL  < 'A-DOLD  ( 2"  (1  ) ,1  =1  ,ME  ) 

I r ( T F AG  , M L • 1 ) WRITE  ( L U N 3 ) ( E D < I > , I = 1 , ME > 
vu;  CONTINUE 

WRIT1'  (ID, 20:2)  C H E 0 ( I > ,1=1*1?) 

N = NFiN  + 1 

w r I t f (I  , ? j ; 7 ) n 

p c T L f i 

2jQ:  format  U6-I0PFMAR7  TAPE  GENERATED  - EXECUTION  TERMINATED* 1 7 H WITH 
i r.Yi;ic  no.  -,ii:,/3.thqthe  current  pestaft  tape  is  tape,i7> 

22'J1  FORMAT  (At  rtQ  RESTART  TAPE  GENERATED  - EXECUTION  CCNT  I CUE  S , 1 7 h W I T h 
10YCL"  NO.  - , 1 1 .1/33  HO  THE  CURRENT  RESTART  TAPE  IS  T ADE  * I 7 ) 

2.12  FORMAT  (IX,  1 2A  <-> ) 

2jD3  F0~--  AT  (4  <♦><  C RE  3 T A -'T  CYCLE  C OMPL  E TE  D- E XE  CUT  lOtl  C OKT  I NUES  , 1 7H  WITH  C 
i y :l:  n . - , i ic  ) 

ENT 


j 


.SIT*1-  JOT  ME  SOL  V'  (P»OP,  0 0?  FF  ,Hl  ,HJ,  Ml  , VJ  , T,  1MI,  JNJ,P,»,  7,C00f , F, 

1 <0,  M , '.MX  ,E0,F°3  ,SIG,  YL  JO  AT,  A ,£,MASS,  F<C  M^R,  TUF  , I FAT  l , 

r ICCM  , 'J 1IM1  , iRTST  , 1 KLOK) 

~ E A l AS  S 

co‘"“*1)M  h -■  s < i ? ) , rri ag, numnp,  numel  ,njmmat  , us  j.*,mt  am , nnn , nt  , t t , crf  t 

1 Np  ? I NT  ,-JP,  Nine.';,  4ST.RA,  Air  a,  JETa,  INTSEC,MflXSEC  ,MTSTOP  , MT3TOP  , I IK 

? mi  , i? , u3 1 ,l  a , I r : , io , luni,  iun2  , i JM3 , lun4 , mm,  lun6 » iun?  , nt  afi 
< m a r, , • i a c i o , lun«  , <r  m 

*’,0'*'*  /l  s a a / l«  ( a) , ss  ( <t , «> , xc,  yc,  s ( o , , c<  4 , a ) , d:  ps  ( A> 

- M/  3 1 *°  J , »K,RL  , ZI  , Z J,  Z<,/L  ,000  ( -*)  ,0T3P  ( «)  , */EL  < «>  , ACCl  ( P ) 
0 I m-'NM  T I PROP  < 7,  1 ) ,COEFF  (6  0 ,1  > , HI  ( i ) , H J (1)  , Jl  ( 1 ) , VJ  (1  ) , T ( 1 > , 

.in:  ( 1 ) , J-J  J ( 1 > , P (?,  1 ) ( 1 ) , Z ( 3 ) , CODE  < i>  , 3(1  ) , xo  ( 1 ) , Xl  ( ‘ ) , X2  ( 1 ) , 

. i x ( ^ , i)  , ps  ( is , i ) , slots , i) , YLraaT ( 3 , i) , a ( i > ,e ( i ) ,mase ( i > , ekc  ( i > , 

.i  ( i ) , -n->  ( i> , oof  ( l ) ,E0  ( l ) 

i I M_ ' S I CN  pc<(  -«) 

DATA  MR  /IM,  hmycle  = ,3,0/ 


i f.1 1 r t at  i x a t ion 

IF  ( I mi  . 0.0)  R'Mf'O  l UNA 
Me  M'UMNMHUMnP 
•1ST  -)l  <=  (HE  . *NE  Q 3-1  ) /ME Qp 
‘ I M 3 -i  L K - ( ’ I U M N p ♦ H 1 - 1 ) / N‘  1 
JEM  iK=  CI'JMEL  ♦ M T-  1 ) /M2 

IF  ( L F L a J . -_Q.  2 ) f:KDl<-  ( MUM:  l fM  J-l)  'M3 

.O'.'MM-,  FOR  Th-  ETTP-OYSTEP  EOLUTIO'I 

A 1 - < . / OH  L T 
A ’-.70/  * L T * *2 
- 2 - A 1 / 2 . 

AS  = ft  2 ♦ A 2 

■ 4 - a i / a j 
as  - ? . / a j 

a (>  - l r / . 

a 7-  L ' * * 2 / 5 . 

A A : A 7 ♦ a 7 

ST  I'TS'ESS  AM  MASS  MATRIX  OF  THE  SYSTEM 

L < \ , 4 * • * 

■ ^ • A j . - J . J 7 

l A l*MANO 

IF  ( I HA IL . - Q. 3 ) SO  TO  4 j 0 
HNM-  ; 

. CCMTINU- 
T C = 'j 
lCpNT  r i 

IF  l I F i AS  . . U.  1 ) SO  TO  5 2 
*EW  I 0 lun 4 
■'EwlNO  MN2 
•’  S W I 0 Q L U i*l  3 
•'E'aIm  ) L* ) N 1 

RE  AC  (1002)  (R(I),I=1»MN) 


,1  E A ° OUMMY  aoc:a 


l JO  SS  I=UNIIMNP 


L ( 2 » T - 1 ) - 3 . 

i ?*  r > = i. 

MASS  ( I ) - 3 . 

•5  3 CONTINUE 

I c < I FL  AG  .Nc  . 1 > GO  TO  b 3 

00  5 7 1 = 1 , LA 

a ( i > - : . : 

5 7 f 1 n n T I n ! I r 
b?  SON-7  T NUC 

^0  780  lMK=  1 ,NEL  JL  K 

1 F ( I r L A G.  NE  . 1 ) F E A C (LUN3  i ( E O ( I ) , I = 1 , M E ) 
NUM/L  0 = 3 

JO  -75  N - 1 , "i 2 
OEL  - M?*  (O'K-1)  +0 
IF  (N'L  . uT  .UU’l-l  ) 5 0 TO  3 7b 
00  1*0  1=1,1* 

J = I +1 

i I = I 3 ( I , • I ) ♦ I X C 1 , N ) 

L 4 ( J ) = I I 

104  LM(J-1)=II-1 


C 

c 


-I*’;  T-*.  CLA3TI.;  CONSTANTS 


i = r x ( i , n ) 

J=I X ( ?,  N> 

< = I v ( 7 , r.  i 
L = I X < •*,  N) 

M T V 1 ) S ( I X ( 5 , N ) ) 

GAL.  C I . >0  (•>,  7 , COOE  , ),  XC,  X 1 , < 2 , I , J,  <,L  ,NCPNT, MNPOLK ,M1 , LUN?  , “X  ) 

:F  ( NON. '!  . . 1 r A.  L STRAIN  ( IX  (1  , N>  ,EOS  ( 1 , N)  , MTV3C  , VOL  , PPOP  (1,  NT  yen  ) 
"ALL  C')N-  CT  ( • PS  ( 1 ,N)  , SI  5(1,0)  , PROP ( 1 , NT  X PC ) , 

Or  EFP  ( l , MT yo  ) , I X ( 1 , N ) ,0,OrLT, NUMVL J, YLOOAT,NEU 


C FORI  ;L  / •''ll  STIFFNESS  MATRICES 

IF ( ‘ 1 0 * 1 . 3. NT)  GO  TO  775 

CALL  SI  IFF ( I F L A i • , L 0 N 4 , M 3 , NUMEL  , NEL , NEQ, VOL  , 

. L X ! 1 , N)  , A^'KC.LK?  ,L  MC,LLMO,PROP  ( 1,  MT  YPE  ) ,M  ASS,  T C,  AO) 

IF  ( NNN . • 0.5)  Ci  0 TO  3 7 5 
30  37.  1=1,4 
3T  7 7 7 J- 1 , 8 
J J = l M l J ) 

7 7:  : ( j j>  = ■ (jj)tss(i,n*siG(i,N)*voL 

S7j  CONTINUE 

3 7c  if  (i  flag,  n . n ^iteiluod  (io(I),i  = i,me» 

4°I'  ( N T A P E ) ( (Cl  >(I,N)  ,1  = 1,4)  , ( E ° S ( I , N > ,1  = 1,4)  » N = 1 , M 2 ) 

if ( n'jmy l j.ve.o)  wpite(Io,2:jo  tt,  ( (ylodat  ( i , j ) , i=i,.  ) , j=i,Nu*n  r) 
2,  .3  FORPAt(?fh  INELASTIC  ELEMENTS  at  T=,  C15.b/8(15H  ElF  FACTOR)/ 

. ( j ( ff  . : , f ■?.  h ) ) ) 

7 8 'j  CONTINC- 

IT  (IFL-.G.rr  .1)  CALL  INTQhG  ( L'J t 1 , L UN 3 ) 

IF  (MNN.EQ.NT)  so  TO  303 

IF ( TFL AG. NC. 1)  CALL  GS T I T F ( L UN4 , L UN 4 , L UN  5 , L UN 1 , NEQB , MO  ANC , N ST P L k , 
.NUMEL  ,M  3,NK>!LK  , A,NEQO,EKC,EKC,  LMC,  LMC,  MASS) 

IF  ( r'NU.NE.  3)  GO  TO  310 
'3  lNITIAi_  ACCELERATION 


I F ( I R l A > . E 0 . 1 ) br  TO  8 30 
REWIMJ  LUN? 


E-  31 


3 E W1 M 0 l JM  * 

i : : jo  <?  ’ : -j  k = i ,hhp  <_  k 

I F (I  F L A G • N .1)  ^FAD  (L(JN2>  <P(I>  ,I  = 1,MN> 

‘IjTanoM'UMMK-lMH 

riST  0°  = ,*'U*N.K 

fjST  ^ p = * I\  0 (NEC,  MS  TOP) 

I K - 1 

*T=r ( 1 , 1) 

"in  h ' r»  i - i , f i 

• • ( 2 * M i > = ; . 

* ( 2 * m 1 - 1 ) - ' . 

z 5 COM  TIN  I £ 

CALL  L )A)(T,P,-‘fr.Jl,JNJ,Hl,HJ,i/Ili/J,I<,NSTAKT,f)STOP,NUMPC,TT,C';LT) 
I I < = I 

jo  f*  i ■;  i - 1 , * i 
jM-j  = •■■!*  cn<-D  * : 

IF (NNF. GT.NUMNP)  jU  TO  515 

v 2 ( ? 7 T - 1 > = A;M1(E»I-1>-E(  2*NNP-1 ) ) / M ASS  ( NN  P > - AL  F A*  X 1 ( 2 * I - 1 ) 

/?(?*!)  - A "j  * ( )(?*!)  - r (2*  0Np  ) >/MASS<NNP>-ALFA*Yl  ( 2*1  ) 
hi,  ;ONT I M" - 

HI'.  IF  ( IF  L V". . 0 . 1)  W-IT£(LUM4)  ( F ( I ) , I = 1 , MN > 

«PTT-  (MAP  .MM  i ) , 1=1,  MN) 
h 2 ;no  T T NO  E 
1 K-I  IK 

IF  ( I F L A , , OF  . 1 ) ALL  TNTChG  (LUN2.LUN4) 

i i . TF  ( JT  . j.  . ) ,0  TO  MS 

CALL  K_C<  ( IKL  n<  , I F TSEC  , MAXSEC  , I )A  I L ,SrC  , I ^EC) 

IF  (OF  ■I.NE.mtsto”)  G:j  T ) 395 
NTS!  TP=  M'-jT  OPFOTi  rC° 

n a i l - 1 

' y IF  ( I 9 A I L • _ G • J ) 0 TO  ♦ 9 5 

M.  CALL  ‘TAT  ( I 0 A I l , LOOM  , P30  P , 4 1 1 '1 1 , * ASS  , c , 0 , MNP  PL  K , - 0 , f 5 L 'L  K , 

1 A , H . T 1 L <) 

I 0 A J '_  = ' 


1 


4 

•4  1 „ 


. j F J = J \ 

: * 1 

I F ( r.  If! 

’ T ) pc  j j p (4 

T T=TT  ♦ 

•EL  r 

IF  (F-mj 

( O'"! 

, i : ) . . . . > 

CALL  MSPLAIPCK 

,NMM) 

1 F ( I r L 

‘ y • z. 

J.l)  Gi  TO 

■♦1C 

•'EM  0- 

L'JN 

7 

00  4 0 9 

n ■>< 

- 1 , f r ml  k 

IF (IF  L 

A 3.N 

1 .1)  PEAL  ( 

lJNZ)  (R ( I ) ,1-1 

, MN  ) 

,'JSTA  * T 

2*  ( 

0 1* ( MOK- 1) ) + 1 

fiSTF  ■»- 

. *‘U 

* N3K 

N S T 0 p - 

M T ' 1 ■ 

1 1 V— 

f NST  CiVJEQ) 

00  4 7 J 

1=1 

* Ml 

' ( j m - 1 ) = : . 
k:*i)  = :. 

OO'lT  t ho  f 
IK- I In 

I AL  L LOAJ  (7  ,P,  9, 1U  1 , JNJ  , Hi  , HJ,  VI  ,\/ J , IK,  NST APT, NSTOP, HUMPC  , T T , CEL  T) 
10  *3C  1=1,  Mi 

MNP=  U*  (tn<-l)  +1 

IF  (HHP.  ,HU‘3N0>  GO  TO  45C 

(,*':np-1)=P(2#I-1)-c.(2*NMH-1)  + M ASS (NNP ) M 14*<1(?*I-1>  ♦ A5*X2(2M-1 
1)  ) 

F.-32 


1 


- ( o*i]  up  ) r 9 (2*1  > - : < 2*  MNP)  + N A SS  (NMP)  * ( A4  * Xi  (2*I>  fAS*X2<  2* I)  ) 

-.3.  CONTINUE 

con r r no  - 
•♦•b ; cnuinuE 
IK-  t < 

If  ( I p L •"  > • E ) • 1 ) ,0  TO  *00 

•-  aI  I i ) t_UN4 

0 0 ‘.-.0  T-1,NST'*LK 

I S T A •>  T - N * H * ( I - : > ► 1 
NS  T 0Q  = N 1 3 * T 

(i  (j)  ,j  = t!STAPT  » N S T 0 P ) 

-j?  0 0'i  T I NO  E 
■♦-jO  CONTI  N j - 

C SAVE  "ACS  VECTOR  FOX  VANEL 

.Ew'TCC  LiJOX 

W p I T ~ ( L IMS)  (MAjSII),I  = 1,0IJMNP) 

CALI  UGCL  1 ( A , P.UF,  MeR,NEQP*  N 3 AN  0 , NST  BL  < , L UN  1 , LUN 5 , L UN4 , L UN6  , LL'K?  , 
. LIJN4  ) 

IF  ( 1FL  TG.  EO.  i>  CO  TO  4b9 
REWIND  LUN4 

jO  4,7  I - l , US  T nL  “■ 

fST7.PT  c J m * C I - 3 > ♦ 1 
f"jTO'J=NE'J-l*  I 

- r ( L : J I 4 ) ( r ( J,,J=j3T7»PT,fvjSTOP) 

4 b 7 CONTI  NOE 

4 •,  ) CON  T T (40  E 
►j op  = o 

I F ( I F L A 0 . EO . 1 ) G TO  4 7 0 

■»  e w r i o l u i j 2 

J F W I NO  l UN 4 

w:  00  49G  N9<=1,HNP-9LK 

1 F f J c L A C , N 1 . 1 ) PEAL  (LUO  2)  (MI)  , I - 1 > ON  > 

N S T 0 p = '**-*  1 

)o  4 / ? i - 1 , nst o •- 

OOr-OOF + 1 

I r ( r of. OT.UFO)  GO  TO  47b 
A C C r A 2 * £ ( N 1 F ) - A I * X 1 (I) - . b * y ? ( I ) 

U I ) = DEL*  *X  i ( 1 ) +A7*  ( ACC  + X2  ( I)  *X2  (I  ) ) 

XO(I)  = nr  ( I ) 4>3(  I > 

*1<I)  - x : ( I ) ♦ A f * ( y 7 ( 1 ) ♦ A cc ) 

/2 ( T ) = ACC 
4 7S  CONTINUE 

47b  ! F ( IFt  A ,.o:  . 1)  rtrvlTF  ( LUN4)  ( R < I ) , I = l , M N ) 

WRIT  I ITA  ->  ■;  ) (X  ( 1 ) , 7 ( i)  , XO  (7*  I - 1 ) ,X0  < 2*1 > * X 1 ( 2*1  - 1)  , X It  2*  I > , 
.'3(2*1-:)  , x ? ( 2 » ! ) , I = 1 , M l ) 

4/)G  C3MTTNUr 

IF  ( I f LA  G.  NE  . i)  '.ALL  iNTCHG  (LON2,lON'+> 

iF(fOT.  j.:>  go  to  i:jo 
I F ( I F L a . . E 0 . 1 ) N 1 = 1 
IF  ( If  L AG.  -E  1.2)  Nl  = KO^f  1 
N2-r  1 + 0 ST 
l3-(  - +NST 
04- N ; *0 ST 
N5=N4  *’  ST  * ? 

Nb  = r. s tN  S’-*  ? 

N 7 -Nb  ♦ 51  * 2 
w«- N 7 40  ST  * ? 

CAlL  VAN^L  ( IXjEC*3.  »Z*COOE*?»XO»Xl»X2»NNP|3L<»NELOL<» 
lAC.il  , A ( N 2 > , A ( N ' ) * A ( N 4 ) , A < NO)  , A ( Nb ) * A ( N 7 ) , UNO  ) ) 

GO  TO  I ; E j 
END 


F-  ( < 


o o 


S-J  T=0..t  i'.r  : T I F P ( 1 FL AG , IT  APT  , IUMEL  , N , k’E  U , V OL  , I x , A , i <C  , L E K , 

: l mt  , Lc ",  m , - as  , t r. , a„> 

C 0 1 ' ■'  1 1 / ^ S 4 A Pi  > / L M ( M ) , S S ( 4 t H ) , X C » V ( ; , 0 ( 3 , 3 ) , ",  ( m , 4 ) , J • ''  S ( 4 ) 

1 i'i  '’I  / f , i 'E  / Ij^Jj^Kj^LjZIiZJj/Kj/LjyCl^fOj  i T S ■*  ( 3 ) 

P-.U  M-;  SS 

TMt  - j . - IV  (l)  ,A(Nr-,D  ,EKC(1)  ,L-I"(i)  ,H(1)  ,Mi3S(l) 

I = I y f 1 ) 

j = I x f ? ) 

< = I X ( 3 ) 

L = I X ( 4 ) 

: r t j . ’ ■ . «)  sc  ’ i:; 

. A l l i'  J f *■  I , J,ZJ,ZJ,3,I/0L) 

',ri  T ° u: 

' A _ L . . AO  ( I , - I , P r , r’L  ,ZI,ZJ,ZK,/L,X'.,Xf  » Z A L ) , c * S > ’ 

i ' ) i r y K'  <*  ze-o  jisplacimemg 

1 A : ■ 1=1,- 

I J-  : * I 

r f-  ( (’)•',•  • J ) >0  T 1 ?UU 

I r‘  ( n r ■ - ( I ) . . . 1 . o ) ,0  T 0 1 b 0 


a < i j , j ) - : . 

; ( i , i j)  - : . 

: • ( sot  e ( i ) . Eu . *: . o ) io  T c ... 

lr  13,  J-1,3 

S ( 1 J - 1 , J>  = . . 


• w 

( J , T J- 

i > = ■: . 

? ] : 

; 0 I T T V ■ j 

z 

r 

r' 

••  L 

E GlOmAL 

MASS  ZfCTQP 

c 

' , A L L r,-i 

S)  i ( I r , M 

ASS, AO , H ( 3) 

c 

C 

AES  < Ml 

; GLO'iAL 

STIFF  JESS  M 

i L 

IF  ( 1 r u 

A . IQ.  2 ) 

>0  TO  uQC 

jO  ' jC 

1-1,3 

II=L  ■*€!  ) 

>o  ,v;o 

1 - « H 

J a.  y 

l 

J J - 1 '.U)-ITM 

IF  ( J I . 

LT.l)  r.i) 

TO  3aT 

A(TI,JJ)=A<JI,JJ 

) ♦ S ( I , J ) 

. j j 

SCOT  I NM 

T 

IF  ( J . L 

T . ( 4 U M c L 1 

T'J^N 

in  ?2: 

1-1, FT 

J- A ( 1 , 1 

) 

if  n . 

J . I . ) GO 

TO  3Z0 

j-  c i * i ) 

/ ’ 
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